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I .  INTRODUCTION 


In  the  last  few  years,  atomic  beam  scattering  has  undergone  a 
period  of  rapid  development.  Beams  of  low  energy  neutral  atoms  are 
used  to  study  the  properties  of  surfaces  a.nd  the  details  of  the  atom- 
surface  i.nteraction .  Although  it  is  possible  to  study  a  variety  of 
phenoir.e.na  such  as  gas-surface  reactions  and  surface  phonon  .modes,  the 
discussion  here  will  be  mainly  restricted  to  elastic  scattering  of 
rare  gas  atoms.  T.he  relative  simplicity  of  this  case  .makes  it  par¬ 
ticularly  useful  for  the  developm.ent  of  exper ime.ntal  and  theoretical 
tec.hniques.  .Moreover,  much  intrinsically  interesting  physics  is 
i.nvolved  in  the  'understanding  of  such  systems. 

The  inform.ation  generated  by  the  atomic  beam  scattering  experi¬ 
ment  is  unique  in  several  ways.  The  lew  energy  and  relatively  large 
mass  of  the  incident  atoms  ass-ure  virfually  no  penetration  into  the 
bulk,  m.aking  the  experiment  alm.ost  totally  S'urface  sensitive.  T.he 
low  energy  of  the  atom.s  (~10  m.eV  or  less  in  some  cases)  makes  them 
an  excellent  probe  of  the  physisorption  regime.  Also,  since  the 
incident  fl'ux  is  quite  low,  it  is  possible  to  gain  direc..  information 
about  the  interaction  of  a  single  atom  in  the  presence  of  a  solid  sur¬ 
face,  as  opposed  to  measurements  of  thermodynam.ic  variables  of 
physisorbed  films  in  which  atom-atom  interactions  must  be  taken  into 
account . 

There  are  also,  however,  several  i.mpcrtant  limitations  in  the 
exper ime.nt.  The  nonzero  velocity  spread  of  t.he  incident  beam  and 
ncnperfect  condition  of  the  sample  surface  tend  to  limit  the 
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resolution.  There  are  no  limitations  in  principle  on  the  type  of 
sample  surface  used,  but  in  practice,  only  a  few  surfaces  have  been  used 
with  much  success.  .Metals,  for  example,  have  produced  diffraction  in 
only  a  few  rare  cases.  Finally,  full  analysis  of  the  scattering  data 
with  physically  realistic  theories  is  extremely  difficult,  especially 
in  the  case  of  inelastic  scattering.  .Much  has  been  done  recently  to 
alleviate  these  problems,  but  they  all  remain, at  least  to  some  extent. 

As  mentioned  above,  the  discussion  here  will  be  mainly  restricted 
to  elastically  scattered  rare  gas  atoms.  Some  of  that  work  has  been 
done  at  relatively  high  energies  (~100  meV) ,  which  primarily  yield 
an  estimate  of  the  corrugation  for  the  periodic  part  of  the  atom- 
surface  potential  by  using  semi -classical  scattering  theory.  This 
study,  however,  will  be  concerned  mostly  with  the  low  energy  regime 
(-10  meV)  where  quantum  effects  predominate.  Under  these  conditions, 
many  proble.ms  of  the  analysis  can  be  avoided  by  considering  only  the 
kinematics  of  the  scattering  event.  This  method  yields  direct  informa¬ 
tion  on  both  bhe  form  of  the  attractive  well  and  the  strongly  repulsive 
periodic  part  of  the  atom-surface  potential. 

The  major  emphasis  of  the  present  wor.k  is  on  the  interaction  of 
helium  atoms  with  the  basal  plane  of  graphite.  An  important  motiva¬ 
tion  for  the  study  of  this  system  is  the  existence  of  much  data  on 
the  thermodynamics  of  submonolayer  He  films  adsorbed  on  graphite. 

These  have  been  widely  used  to  e.xplore  the  physics  of  two-dimensional 
arrays,  a.nd  it  is  hoped  that  the  present  work  will  contribute  to  t.he 
understanding  of  this  body  of  work. 
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II .  HISTORICAL  REVIEW 

A.  Early  Work 

The  earliest  work  on  atomic  beam  scattering  from  solid  surfaces 
is  that  of  Stem,  et  al?"  in  1930.  They  observed  diffraction  of  He 
and  H^  atoms  from  alkali  halide  surfaces,  thus  confirming  the  de  Broglie 
wave  hypothesis.  They  also  saw  anomalous  mini.ma  in  the  specular  beam 
intensity.  Lennard-Jones  and  Devonshire^  hypothesized  that  these 
minima  were  evidence  of  the  atoms  entering  into  bound  states  of  an 
attractive  well  between  the  atom  and  surface.  They  termed  this  phe¬ 
nomenon  "selective  adsorption". 

The  beam  used  in  this  wor.k  was  an  "oven  beam"  with  a  Maxwell- 

Boltzman  distribution  of  velocities.  This  gave  poor  resolution.  The 

vacuum  attainable  at  that  time  was  also  poor.  Rapid  advances  in 

vacuum  technology  led  to  a  renewal  of  interest  in  beam  scattering  from 

solids  in  the  mid  1960 's.  Owing  to  the  continued  use  of  thermal  beams, 

3  4 

however,  the  resolution  obtained  by  Crews,  Fisher  and  Bledsoe,  and 
O'Keefe  et  al  was  little  better  than  Stern's.  These  workers  all 
used  alkali  halide  targets;  O'Keefe  et  al .  tried  different  isotopic 
forms  (^He,  D^)  and  heavy  rare  gases  for  incident  beams. 

B .  Recent  Work 

The  next  major  advance  was  the  use  of  "supersonic  nozzles"  to 

.  6 

obtain  very  narrow  velocity  distributions.  Using  nozzle  beans  and 

ultra-high  vacuum  techniques,  accurate  studies  of  ‘*He  diffraction  from 

7  8 

alkali  halides  were  done  by  Williams  et  al.  and  Boato  et  al.  H  and 
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D  were  scattered  from  alkali  halides  in  ultra-high  vacuum  by  Wilsch 
9 

et  al.,  but  oven  beams  had  to  be  used  to  keep  the  atomic  species 

dissociated.  These  workers  all  cleaned  the  surface  by  heating  in 

10 

vacuum.  Houston,  .Meyers,  and  Frankl  obtained  a  clean  surface  by 
cleaving  the  alkali  halide  crystals  ^  situ  and  studying  '*He  scatter¬ 
ing  from  them.  Williams  et  al.  and  3oato  et  al .  concentrated  mostly 
on  diffraction  work, while  Wilsch  et  al.  and  Meyers  et  al.  emphasized 
study  of  the  bound-state  resonances  (selective  adsorption) . 

With  the  advent  of  these  higher  quality  data,  new  theoretical 

efforts  were  undertaken.  The  early  work  of  Cabrera  et  al.^^  and  of 
12 

Wolken  did  not  properly  include  the  resonance  states  and  predicted 
maxima  in  the  specular  beam  instead  of  the  experimentally  observed 
minima.  The  calculations  of  Chow  and  Thompson^^  resolved  this  problem, 
predicting  specular  minima  within  an  elastic  scattering  framework. 
Another  problem  with  older  theories  was  the  use  of  unrealistic  model 
potentials,  e.g.,the  .Morse  potential,  which  was  shown  to  be  incon¬ 
sistent  with  the  more  accurate  data.^^  Pairwise  sums  and  semi- 

14 

empirical  potentials  have  since  been  employed,  along  with  simplified 
models  {e.g.,a  hard  wall  with  corrugation  and  square  well)  in  some  of 
the  scattering  calculations. 

The  most  recent  theoretical  work  has  attempted  to  account  for  the 

fine  structure  seen  in  atom-surface  scattering.  Examples  of  this  are 

the  work  of  Weare  et  al.^^  and  Garcia  et  al.^^  Theoretical  accounts 

17  18 

of  the  interactions  of  bound  state  resonances  '  have  also  been 
attempted.  The  latter  phenomenon  has  been  t.he  subject  of  recent 


experimental  investigation  also.  Another  trend  in  recent  experi¬ 
mental  work  has  been  toward  a  wider  diversity  of  surfaces  used  as 

20  I* 

targets.  Cardillo  has  seen  diffraction  of  He  from  the  (111)  sur- 

21  22 

face  of  silicon.  Work  has  also  been  reported  on  silver,  nickel, 

and  nic.kel  oxide,  using  beams  of  ‘*He  and  H2  •  T.he  basal  plane  sur- 

24 

face  of  graphite  has  been  extensively  studied  by  Boato  et  al.,  by 

25  26 

Valbusa  et  al.  using  Hj »  ^2'  author.  Boato 

et  al.,  using  beams  of  ‘‘He  atoms  with  energies  of  -17  meV  and  ~63  meV, 


made  a  thorough  study  of  the  interactions  between  resonances  and  a 
preliminary  measurement  of  the  resonance  binding  energies,  as  well 
as  some  diffraction  work.  The  contributions  of  the  author  are 
discussed  in  the  remainder  of  the  present  thesis. 
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III .  THEORY 

A.  Introduction 

The  development  of  a  complete  and  satisfactory  theory  for  atomic 
beam  scattering  is  far  from  finished.  Great  strides  have  been  made 
in  recent  years,  however.  More  importantly,  a  large  class  of  experi¬ 
ments  can  be  interpreted  extremely  well  using  a  set  of  simplifying 
assumptions.  The  latter  will  first  be  developed  extensively,  and 
then  some  attention  will  be  devoted  to  the  former. 

A  complete  scattering  theory  must,  of  course,  take  into  account 
the  presence  of  inelastic  events  in  the  scattering  process.  At  pre¬ 
sent,  however,  such  a  theory  is  virtually  nonexistent  and,  moreover, 
has  little  relevance  to  most  of  the  experimental  work  described  here¬ 
in.  Therefore,  we  will  restrict  the  discussion  to  elastic  scattering. 
The  general  goal  of  these  theories  is  to  predict  the  intensities  of 
the  scattered  beams.  It  is  possible,  however,  to  obtain  much  useful 
information  solely  from  the  kinematics  of  the  scattering  event  with¬ 
out  regard  for  the  actual  intensities  of  the  outgoing  beams;  only 
their  geometry  with  respect  to  the  surface,  along  with  that  of  the 
incident  beam,  need  be  known. 

Only  a  brief  summary  of  theoretical  considerations  is  given  here 
since  the  author  has  not  made  a  significant  contribution  in  this 
area.  The  main  purpose  of  this  chapter  is  to  clarify  the  interpreta¬ 
tion  of  the  data.  In  Sections  B  and  C, the  kinematical  point  of  view, 

2 

first  developed  by  Lennard-Jones  and  Devonshire,  is  discussed.  In 
Section  D, a  brief  examination  of  several  recent  dynamical  theories  is 


presented. 
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B.  Kinematics  of  the  Scattering  Process 

The  incident  atom  is  well  modelled  as  a  free  particle  character¬ 
ized  by  a  kinetic  energy,  This  is  essentially  a  plane  wave  of 

wave  vector  j  |  =  ^  ,  and  the  incident  wave  vector  can  be  decomposed 
into  components  parallel  and  perpendicular  to  the  sample  surface: 


=  ic  -t-  k^  2  (la) 

and 

K  =  k  X  +  k  y  .  (lb) 

X  y 


The  z-axis,  of  course,  is  the  surface  normal,  while  the  x  can  be  chosen 
to  be  parallel  to  some  crystal  symmetry  axis  in  the  surface  lattice. 


The  two-dimensional  reciprocal  lattice  of  the  crystal  surface  can 
be  obtained  via  the  usual  techniques  from  the  corresponding  two- 
dimensional  real-space  lattice  of  the  surface.  Vectors  in  this  space 
are  denoted  as 

=m,n  =  ^  .  (2) 

where  g^  and  g  are  the  fundamental  vectors  of  the  surface  reciprocal 
lattice . 


In  order  to  satisfy  conservation  of  crystal  momentum  in  the 
scattering  event,  the  two-dimensional  Bragg  condition 


K 


f 


K. 

1 


+  G 

m,n 


(3) 


must  be  satisfied  for  the  final  and  initial  wave  vector  components 
parallel  to  the  surface.  Conservation  of  energy  can  be  ensured  by 
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requiring  that 


^  .  2  _  „ 
2in  i  ^i 


=  E(K^)  +  E. 


zr 


(4) 


for  the  final  and  incident.  Note  that  these  conditions  are  only  valia 
for  the  case  of  elastic  scattering. 

Combining  Equations  G) and  (4) and  approximating  the  atoms  trans¬ 
lational  motion  as  a  free  particle  yields  a  condition  for  an  elastic 
scattering  event  to  occur;  namely, 


.  +  G 
1  m,n 


+  k 


2 

Zf 


(5) 


If  k^^  is  positive, there  will  be  an  outgoing  scattered  beam;  this  is 
simply  the  case  of  diffraction  from  a  two-dimensional  lattice.  If 
k^^  is  negative, there  will  be  no  diffracted  beams  allowed.  There  is 
a  physical  situation,  however,  which  corresponds  to  negative  k^^. 

Van  der  Waals  forces  give  rise  to  a  weak  attractive  potential 

between  the  gas  atom  and  the  surface,  leading  to  a  discrete  set  of 

2 

bound  states  with  energy  E.  =  -r—  C.  .  when  k  .  is  equal  to  one  of 

3  2m  3  zt 

the  €. ,  the  incident  atom  is  in  resonance  with  a  bound  state  and  under- 


27 

goes  a  so-called  "selective  adsorption"  transition.  The  selec¬ 

tively  adsorbed  atom  is  in  a  one-dimensional  bound  state  with  respect 
to  the  z-axis  and  translates  across  the  surface  in  a  Bloch  state. 


Its  translational  motion  can  be  approximated  as  a  free  particle,  in 

•ft  2 

which  case  its  translational  energy  is  simply  q—  K,.  In  this 

2m  E 


approximation,  then,  the  selective  adsorption  condition  becomes 


(6) 


(K,  > 


in,n 


=  k?  - 


t.  . 

3 


= 


This  IS  readily  seen  to  be  the  equation  of  a  circle  in  space 

centered  at  -G  and  having  a  radius  of  (k?  + 

m,n  ^  1  '  3  ' 

Notice  that  the  are  the  eigenvalues  of  a  one-dimensional 
Schroedinger  equation  employing  the  laterally  averaged  atom-surface 
potential  Hamiltonian.  Thus,  much  information  about  the  nature 
of  this  potential  can  be  extracted  from  measurements  of  the  .  These 
in  turn  can  be  calculated  from  Equation  (6)  if  the  magnitude  of  the 
incident  wave  vector  and  its  orientation  relative  to  the  surface  are 
known  when  a  selective  adsorption  occurs.  The  effects  of  selective 
adsorption  on  the  scattered  intensity  will  be  discassed  in  Section 
II. D.  The  reliability  of  the  approximation  leading  to  Equation  (6) 
will  be  assessed  in  the  following  section. 


C.  Band  Structure  of  Adsorbed  Atoms 

The  potential  energy  of  an  atom  in  the  presence  of  a  solid  surface 
can  be  written  as  a  Fourier  expansion  in  the  surface  reciprocal  lattice 
vectors: 

00  *  . 

V(r)  =  Vg(z)  +  ^  V-  (z)e^^’^  .  (7) 

The  (z)  term  is  the  laterally  averaged  potential  and  the  summation 
of  higher-order  terms  constitutes  the  two-dimensional  periodicity  of 
the  surface,  i.e.,the  corr-jgation.  The  free-atom  approximation  made 
in  the  preceding  section  consists  essentially  of  retaining  only  the 


G  =  0  term  of  the  expansion.  The  approximation  may  be  refined  by 
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including  the  higher-order  terns  as  a  perturbation  or  by  including 
them  in  dynamical  scattering  calculation.  Using  the  latter  technique, 
Chow  and  Thompson^ ^  predicted  the  effects  of  higher-order  terms.  If 
circles  predicted  by  Equation  (6)  for  two  different  states  [  labelled  by 
Cm,n)  and  j]  cross  in  space, there  will  be  a  degeneracy  of  the 

two  states  there.  The  periodic  part  of  the  potential  acts  to  lift  the 
degeneracy,  giving  rise  to  a  band  splitting.  This  is  observed  in  the 
scattering  data  as  a  splitting  of  the  selective  adsorption  features 
instead  of  their  merging  into  a  single  feature. 

The  magnitude  of  the  splitting  depends,  of  course,  on  the  Fourier 

component  which  admixes  the  states  as  well  as  their  energy  quantum 

7  8 

numbers.  Standard  second-order  perturbation  theory  yields  the 

result  that  the  magnitude  of  the  splitting  is  equal  to  2  <jiv- 

t.he  matrix  element  which  connects  the  two  states.  In  many  cases,  this 

matrix  element  is  too  small  to  produce  an  observable  splitting  and  the 

levels  appear  to  cross.  A  more  sophisticated  band-structure  calcula- 

29 

tion  has  been  done  by  Carlos  for  He/graphite, and  it  was  found  that 
the  perturbation  theory  results  are  generally  quite  good.  Although 
there  are  still  small  deviations  in  regions  far  from  crossings,  the 
behavior  there  is  quite  close  to  free-atom-like. 

The  foregoing  discussion  has  several  implications  for  the 
experiments.  In  order  to  accurately  measure  the  bi.nding  energies 
of  the  physisorbed  atoms,  the  data  must  be  taken  far  from  any  split¬ 
tings  ,  which  render  Equation  (6)  incorrect.  If,  however,  data  is  taken 


extensively  throughout  a  splitting,  the  matrix  element  for  that  case 


11 


can  be  calculated.  Thus,  information  can  be  obtained  concerning  the 
periodic  part  of  the  potential  as  well  as  (z)  ,  which  can  be  recon¬ 
structed  from  the  E . . 
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Dvnamical  Scattering  Theories 


There  are  many  problems  inherent  in  the  prediction  of  actual 
intensities  of  scattered  beams  as  a  function  of  incident  geometry  and 
energy.  Some  of  these  are  experimental,  such  as  the  imperfection  of 
the  surface  and  distribution  of  energies  in  the  incident  be--?.m.  A 
major  theoretical  problem  is  the  inability  to  handle  inelastic 
scattering,  which  is  often  quite  important.  Other  problems  are  the 
severity  of  the  approximations  necessary  to  make  the  calculations 
tractable  and  nonuniqueness  of  the  models  which  provide  good  fits 
to  the  data. 


Despite  these  problems,  much  progress  has  been  made  in  recent 
years.  The  semiclassical  approaches  will  not  be  dealt  with  here 
since  they  are  not  very  relevant  to  the  energy  regime  of  interest.  A 
good  review  of  these  methods  is  given  by  Goodman. Two  of  the  more 
fruitful  techniques  undertaken  recently  will  be  discussed. 

The  method  of  Chow  and  Thompson^^ ' is  to  expand  the  atomic  wave 
function  in  surface  reciprocal  lattice  vectors  and  insert  it  into  the 
Schroedinger  equation  along  with  Equation  (7).  This  yields  a  set  of 
coupled  differential  equations: 


i 


which  are  numerically  integrated.  Of  this  infinite  set  of  equations, 

some  set  N  must  be  chosen  for  the  computation,  corresponding  to  the 

N  diffraction  channels  (i.a.,G  vectors)  which  are  important  in  the 

scattering  process.  It  was  discovered  that  certain  closed  channels 

(k",<0)  must  be  included  as  well  as  ooen  channels.  When  this  was  done 
zf 

properly, the  theory  predicted  intensity  mini.ma  for  the  selective 
adsorption  features  in  the  specular  beam,  in  agreement  with  the 
experimental  results.  This  was  found  to  be  the  case  for  several 
different  model  potentials. 

A  Green's  function  formalism  was  m.ore  recently  employed  by  Weare 
et  al.^^  to  predict  the  lineshapes  of  the  resonance  features.  They 
write  the  scattering  equation  in  integral  form  : 


where  Ixn^  is  an  eigenfunction  of  the  unperturbed  Hamiltonian, 
the  periodic  part  of  the  potential  [  given  by  the  summation  in 


(9) 


V  is 
P 

Equation 


(7)  ],and  is  a  representation  of  Green's  operator.  This  method  leads 
to  singularities  at  the  resonance  energies.  These  are  eliminated  by 
means  of  projection  operators.  Approximate  expansions  and  numerical 
techniques  are  then  used  to  solve  the  scattering  problem  in  detail. 


The  major  qualitative  results  of  this  method  are  readily  summa¬ 
rized.  Concerning  the  lineshape  of  the  selective  adsorption  feature, 
a  minimum  in  specular  intensity  is  expected  if  the  resonant  state 
is  strongly  coupled  to  both  the  incident  beam  and  an  open  diffraction 


channel.  Othenvise,  a  ma.ximum  in  intensity  is  expected  to  be  the 
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dominant  feature.  Another  interesting  prediction  concerns  resonances 
which  are  not  observable  due  to  weak  coupling  or  high  energy  quantum 
number.  Such  a  state  could  become  observable  by  "borrowing"  intensity 
from  a  more  prominent  resonance  to  which  it  is  strongly  coupled  by  a 
low-order  Fourier  component.  This  would  happen  near  crossi.ngs,  as 
discussed  in  Section  III.C. 

A  set  of  rules  to  predict  the  lineshape  of  the  specular  intensity 

16 

at  resonance  has  also  been  derived  by  Garcia  et  al .  using  a  corru¬ 
gated  hard-wall  with  an  attractive  well  as  a  model  potential.  These 
calculations  employ  an  exact  scattering  formalism  and  yield  results 
similar  to  those  of  Weare  et  al.,^^  although  the  rules  are  stated  in 
different  terms.  This  .method  can  also  be  used  to  treat  the  mixing  of 
resonances.  The  predictions  of  these  theories  are  compared  with 
experiments  i.n  Chapters  V  and  VI  after  the  experimental  results  have 
been  presented. 
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IV.  APPARATUS 

The  xa]or  conpor.e.nts  of  the  atonio  beam  apparatus  were  constructed 
31  32 

by  D.  E.  Houston  and  J.  E.  Meyers.  A  description  of  the  essential 
parts  of  the  apparatus  is  given  in  Sections -A  and  3.  The  author  also 
made  several  contributions  to  and  modifications  of  the  equipment. 

In  Section  A,  the  investigation  of  beam  instability  and  the  brief  men¬ 
tion  of  low-temperature  .nozzle  cooling  belong  in  this  category.  In 
Section  3,  the  modifications  to  the  sample  holder  are  the  major  example. 

The  work,  described  i.n  Sections  C,  D,  E,  and  F  was  done  i.n  whole 
or  i.n  part  by  the  author.  In  Section  C,  the  problems  of  apparatus 
alignment  are  considered.  Section  D  describes  the  auxiliary  system 
devised  in  order  to  perform  ^He  scattering  experi.ments .  The  equip- 
m.ent  and  teohr.iques  needed  for  cryogenic  cooling  of  the  nozzle  are 
discussed  in  Section  E,  while  Section  F  presents  a  brief  account  of 
the  electron  optics  and  its  uses. 

A.  Beam  Formation 

The  first  step  in  the  form.ation  of  the  atomic  beam  is  e.xpansion  of 
the  source  gas  through  a  .nozzle  into  vacuum.  The  nozzles  used  in  t.his 
work  are  electron  m.icroscope  apertures  (platinum  or  m.olybde.num) , 
typically  about  10  microns  in  diam.eter.  The  second  part  of  th  ’  beam 
source  is  an  electroform.ed  copper  ccr;e  ("skimmer")  with  an  orifice 
of  about  0.2  mm  dia.meter  facing  the  nozzle  as  shown  in  Figure  1.  The 


skimm.er  can  be  moved  laterally  for  alignm.ent  purposes. 
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The  details  of  supersonic  beam  formation  have  been  discussed  at 
6,27 

length  elsewhere,  and  only  the  more  important  properties  of  nozzle 

beams  will  be  summarized  here.  The  motion  of  the  gas  atoms  becomes 
more  directed  along  the  axis  of  the  nozzle  aperture  with  the  skimmer 
eliminating  most  of  the  remaining  gas.  The  gas  temperature  decreases 
greatly  during  the  expansion,  but  the  average  kinetic  energy  of  the 
atoms  comprising  the  beam  is  determined  by  the  temperature  of  the  gas 
behind  the  nozzle.  The  distribution  of  velocities  is  much  narrower 
than  a  .Maxwellian  distribution  characterized  by  the  temperature  of  the 
gas  behind  the  nozzle.  The  details  of  the  velocity  distribution 
depend  on  the  source  pressure,  nozzle  diameter,  and  nozzle  tempera- 
ture . 

In  general,  the  distribution  narrows  with  increasing  pressure 
and  diameter.  With  the  nozzle  cooled  by  liquid  nitrogen,  typical 

'  V 

source  pressures  (-6  atm)  yield  velocity  distributions  of  ^  <2% 

FWiM.  At  a  nozzle  temperature  of  20K,  typical  source  pressures  are 
lower  (~1  atm)  and  the  distribution  is  a  little  broader  —  "2  4  ^ 

2  V  2 

FVJHMj .  The  upper  limit  of  the  source  pressure  is  essentially  deter¬ 
mined  by  the  pumping  speed  in  the  first  stage,  because  if  the  pressure 
there  becomes  too  high, the  beam  intensity  attenuates  rapidly. 

The  skimmer  separates  the  first  stage  from  the  second  stage 
(see  Figure  2).  The  second  stage  houses  the  mec.hanical  chopper. 

The  beam  is  chopped  at  240  Hz,  while  a  light  and  photo  cell  provide 
a  reference  signal  to  a  lock-in  amplifier.  This  allows  e.xtraction 
of  the  bean  signal  from  the  background  noise,  which  may  be  several 
orders  of  magnitude  larger. 


Flexible 
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The  third  stage  is  a  velocity  selector  consisting  of  five  rotating 
32  34 

slotted  disks.  '  This  is  needed  only  for  neasuring  the  velocity 
since  the  velocity  distribution  is  very  narrow.  During  the  scattering 
experiments,  the  selector  is  lifted  out  of  the  beam.  When  measuring 
the  velocity,  the  selector  is  run  in  both  directions  and  an  average 
taken  to  compensate  for  small  misalignments  in  the  rotation  axis. 

The  fourth  stage  is  a  buffer  for  the  ultrahigh  vacuum  in  the 
scattering  chamber.  Separating  them  is  a  collimating  pinhole  which 
is  the  defining  aperture  tor  the  beam  incident  on  the  sample.  The 
long  distance  traveled  by  the  beam  (5:90  cm)  and  the  narrowness  of  the 
pinhole  (fO.3  mm)  result  in  a  highly  collimated  beam  (19  <  10  ’  rad). 

The  major  problem  wit.h  the  beam  formation  system  has  been  the 
occurrence  of  momentary  pressure  bursts,  especially  in  stage  two. 

These  result  in  anomalous  drops  in  the  beam  intensity.  The  problem 
has  been  partially  alleviated  by  installing  a  larger  second  stage 
diffusion  pump,  regulating  the  cooling  water  flows  to  all  the  diffu¬ 
sion  pumps,  and  avoidi.ng  contaminated  0-rings.  The  problem  still 
exists  to  some  extent,  however.  Small  pressure  fluctuations  also 
cause  some  noise  in  the  beam  intensity. 

B .  Scattering  Chamber 

The  beam  formation  system  is  connected  to  the  scattering  chamber 
by  a  flexible  bronze  bellows  to  allow  positioning  of  the  beam.  The 
chamber  is  pumped  by  a  turbo-molecular  pump  backed  by  a  small  diffusion 


pump  to  aid  in  pumping  light  gases.  A  titanium  sublimation  pump  is 


also  used,  although  it  is  not  needed  to  maintain  the  ultimate  pres¬ 
sure  once  it  has  been  realized.  The  system  is  baked  at  ~100°C  and 
this  yields  a  nominal  pressure  of  =3  x  10“^°  torr.  The  titaniiim 
sublimation  pump  is  mounted  in  a  Tee  equipped  with  a  jacket  which  is 
cooled  with  liquid  nitrogen  to  condense  contaminants  whenever  the 
sample  surface  is  cooled.  As  shown  in  Figure  3,  the  scattering  chambe 
has  a  number  of  ports  for  the  beam,  sample  manipulators,  windows,  and 
electron  optics.  These  ports  are  knife-edge  flanges  sealed  by  copper 
gaskets.  A  large  flange  on  top  has  a  rotary  feed-through  for  the 
detector  and  a  port  for  the  ionization  gauge.  The  bottom  flange  has 
a  rotary  feed-through  for  the  sample  holder  motion,  electrical  feed¬ 
throughs,  and  a  sealed  liquid  nitrogen  inlet. 

The  sample  holder  is  constructed  so  as  to  allow  rotation  of  the 
sample  about  two  mutually  perpendicular  axes,  as  shown  in  Figure  4. 

The  holder  is  mounted  on  a  rigid  shaft  which  extends  to  outside  of 
the  chamber.  Three  teflon  rings  divide  the  feed-through  area  into 
two  differentially  pumped  regions.  Thus,  the  polar  angle  can  be 
read  directly  on  a  vernier  scale  mounted  on  the  outside.  Changes  in 
angle  can  be  measured  to  a  resolution  of  0.1°,  but  absolute  angie 
measurements  are  more  difficult  (see  Section  IV. C.).  A  mechanical 
rotary  feed-through  drives  a  chain  and  sprocket  to  rotate  the  bearing 
mounted  inner  section  of  the  sample  holder.  This  arrangement  has  a 
good  deal  more  backlash,  but  the  azimuthal  angle  scale  is  mounted 
directly  on  the  holder  and  read  through  a  window  to  insure  accuracy. 
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Figure  3.  Top  view  of  scattering  chamber. 
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It  is  often  necessary  to  heat  or  cool  the  sample.  Liquid  nitro¬ 
gen  coolant  is  fed  from  outside  through  flexible  bellows-type  copper 
tubing  into  a  copper  reservoir  block.  This  block  is  then  maintained 
in  good  thermal  contact  with  the  sample  holding  mecha.nism.  In  work 
with  the  alkali  halides,  a  rather  massive  sample  holder  is  needed  for 
cleaving  in  vacuum.  Thermal  contact  is  made  in  this  case  by  sliding 
copper  plates.  The  sample  could  be  heated  to  -450  K  by  a  tantalum 
wire  resistive  heater  encased  in  ceramic  tubing. 

For  work  with  graphite  samples,  the  large  mass  is  not  needed,  even 

for  the  vacuum  cleaved  surfaces  (see  Section  V.3.).  Also,  it  is  more 

common  to  heat  the  sample  to  high  temperatures  (~700  K)  to  clean  the 

surface,  a  technique  which  requires  low  mass  and  good  thermal  isolation. 

This  was  accomplished  by  attaching  a  lightweight  copper  sample  mount 

to  a  copper  shaft  which  is  held  in  suspension  by  eight  thin  stainless 
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steel  screws.  Emoedded  in  the  sample  mount  is  a  Spectra-mat  heater. 
Thermal  contact  with  the  liquid  nitrogen  reservoir  is  made  via  a  flex¬ 
ible  copper  braid;  inserted  between  the  braid  and  the  shaft  is  a 
beryllia  tab  which  conducts  well  at  low  temperatures  for  cooling  a.nd 

insulates  at  high  temperatures  for  thermal  isolation.  The  sample  can 
be  cooled  to  -100  K  and  heated  to  -700  K  as  measured  by  a  chromel- 
alumel  thermocouple  mounted  on  the  shaft. 

A  later  refinement  was  the  replacement  of  the  eight  thin  mounting 
screws  by  three  screws  through  a  front,  spring-loaded  rocker  plate. 

This  arrangement,  shown  in  Figure  5,  allows  the  alignment  of  the 
azimuthal  axis  perpendicular  to  the  polar  axis  and  t.he  beam.  Align¬ 
ment  is  discussed  further  in  the  followi.ng  section. 
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The  detector  is  mounted  on  a  rigid  shaft  with  a  feed-through 
similar  to  the  one  on  the  sample  holder.  The  detector  swings  along 
an  arc  with  =50  mm  radius,  and  it  has  a  vertical  travel  of  -5  cm 
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achieved  by  running  the  shaft  through  a  fle.xible  sealed  bellows. 

This  allows  the  detector  height  to  be  adjusted  for  opti.mum  bea.m 
intensity  and  makes  possible  work  with  out-of-plane  diffracted  beams. 
The  angle  scale  of  the  detector  is  mounted  outside  the  vacuum  and 
can  be  read  to  0.1'"  (see  Section  IV. C.);  changes  in  detector  height 
can  be  .measured  to  0.01  mm. 

The  detector  is  a  quadrupole  mass  spectrometer  tuned  to  the  mass 
of  the  incident  atom  (“He  or  ^He)  .  The  bea.m  enters  the  ionizer  box 
through  a  narrow  (~0.a  mm  x  3  mm)  aperture.  The  ionized  He  ator.s  are 
selected  out  by  the  quadrupole  filter  and  detected  by  an  electron 
multiplier.  The  signal  from  the  electron  multiplier  is  fed  via  a 
current-sensitive  preamplifier  to  the  input  of  a  lock-in  amplifier 
provided  with  the  chopper  reference  signal.  The  output  of  the  lock- 
in  amplifier  is  proportional  to  the  intensity  of  the  bean,  but  it  is 
not  calibrated  to  an  actual  flux  since  the  intensities  of  the  scat¬ 
tered  beams  can  be  normalized  to  incident  intensity.  Furtherm.ore , 
many  interesting  phenomena  depend  only  on  relative  changes  in  the 
scattered  intensities. 

C .  Alignment 

Although  the  goniometers  of  the  sample  and  detector  can  be  read 
to  accuracies  of  ~0.l’,  true  angle  readings  of  this  accuracy  can  only 
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be  performed  if  the  apparatus  is  properly  aligned.  Ideally,  the  fol¬ 
lowing  conditions  should  exist:  The  polar  angle  rotation  axis  and  the 
detector  rotation  axis  should  be  collinear.  They  should  also  be 
perpendicular  to  the  beam  and  cross  it  at  some  point.  T.he  azinut.hal 
axis  should  cross  t.hrough  the  same  point  and  be  perpendicular  to  both 
t.he  bea.Ti  and  t.he  other  axes  when  c  =  0°  .  Finally,  also  when  c  =  0'^  , 
the  Surface  should  lie  on  the  plane  defined  by  the  beam  and  the  polar 
axis  with  its  normal  collinear  to  the  azimuthal  axis. 

In  practice,  .none  of  these  conditions  can  be  met  perfectly.  If 
extreme  care  is  taken  in  leveling  and  positioning  t.he  beam  (using 
finely  threaded  screws),  aligning  the  axes  (done  optically , using 
shi.ms  and  the  rocker  place  described  above)  ,  and  positioning  the 
sample,  the  misalign.ments  can  be  reduced  to  about  'l""  in  the  angles 
and  '0.1  mm.  in  the  distances.  The  actual  situation,  however,  is 
worse  than  t,his  implies  because  the  crystal  surfaces  are  net  flat 
planes  but  rather  sets  of  misaligned  facets. 

The  solution  to  this  latter  problem  depends  on  the  nature  of  ti'.e 
sam.ple  surface  bei.ng  studied,  graphite  being  the  mo.st  intractable. 
Therefore,  these  cases  will  be  discussed  individually  in  Section  ”. 

In  each  case,  minimizing  the  mj-salignm.ents  described  here  is  pre¬ 
requisite  to  obtai.ning  accurate  results.  A  detailed  exposition  of 
the  effects  of  misalign.ments  and  tec.hniques  chosen  to  deal  with  tr.em 
IS  given  in  the  .Appendix. 


D.  He  Recirculation  System 

Most  of  the  gas  coring  through  the  nozzle  is  retained  in  the  first 
stage  by  the  skiniiner  and  pumped  away.  In  the  case  of  '*Ke,  the  source 
gas  is  provided  by  a  high-pressure  cylinder  of  research  grade  (99.9995% 
pure)  helium  and  the  waste  gas  is  simply  vented  to  the  atmosphere 
through  the  first  stage  forepump.  This  method  cannot  be  employed  for 
^He  due  to  its  high  cost.  Instead,  a  system  was  constructed  to  recir¬ 
culate  the  used  gas  (except  the  small  amount  which  goes  through  the 
sk  1. timer )  .  T.he  gas  is  pu.m.ped  into  a  reser'/oir  for  storage  v'  in  not  in 
use , 

A  schematic  representation  of  this  syste.m  is  sho'/m  in  Figure  6. 

The  vent  of  the  forepump  is  sealed  from  atmosphere  with  a  valve, and 

the  exhaust  gases  are  fed  through  a  zeolite  trap  immersed  in  licuid 
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nitrogen.  The  'He  is  t.nen  further  purified  by  a  Randex  inert  gas 

purifier  ccnsisting  of  a  titanium,  "sponge"  heated  to  1100  K  inside 
a  guartz  tube.  The  cure  gas  is  then  pressurized  by  a  sealed  refrigera¬ 
tor  co.mpressor  capable  of  producing  a  source  pressure  of  -8  atm..  The 
earlier  alkali  halide  work  was  done  with  a  sealed  diac.nragm^^  ccm.cressc; 
t.hat  only  produced  at.m  with  the  am.ount  of  ^He  available.  The 
pressurized  gas,  of  course,  is  fed  to  the  nozzle  ccm.pleting  the  cycle. 

This  recirculation  system,  is  interfaced  with  the  ^He  source  and 
the  ''rie  storage  cylinders  by  system.s  or  valves  and  pum.ps.  '.'I'lth  the 
valve  separating  the  ccm.pressor  from,  the  nozzle  closed,  the  **'rie  is 
evacuated  from,  the  volume  behind  the  nozzle.  The  valve  to  the  "He 
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source  is  then  closed  and  the  valve  to  high-pressure  output  of  the  com¬ 
pressor  opened.  The  pressure  is  regulated  by  a  precisio.n  feedback 
valve  between  the  compressor's  output  and  input  sides.  When  the  run 
is  completed,  this  valve  is  opened  completely  and  the  ^rie  is  pumped 
out  by  a  sealed  .mec.hanical  pump  with  its  vent  connected  to  t.he  storage 
cylinders.  3y  manipulating  the  valves  shown  in  Figure  6,  the  same 
pump  is  used  to  transfer  the  ^He  from  the  storage  volume  into  the 
recirculation  system.  Fi.nally,  a  mercury-filled  glass  Toepler  pump 
is  used  to  tra.nsfer  the  ^He  from,  the  rather  large  vol'um.e  at  the  ve.nt 
of  the  .mecha.nical  pu.mp  to  either  the  storage  cylinders  or  recirculation 
syste.m. 

Great  care  must  be  taken  to  insure  that  no  lea.kage  of  He  occurs. 

All  devices  are  sealed  and  lea.k-tested,  the  valves  are  bellows  sealed, 
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and  pl’umDing  connections  are  either  brazed  or  swageloc.ked .  The  loss 
rate  of  the  ^Ke  is  knowTi  empirically  to  be  '■0.031  l-at.m/hr,  but  it 
is  not  known  how  much  of  t.nis  gees  through  t.he  ski.mm.er.  The  m.axi.mum. 
pressure  generated  by  the  ccm.presscr  is  a  limiting  factor  in  achieving 
a  narrow  velocity  distribution.  This  is  particularly  a  prcble.m  'with 
He  si.nce  it  has  a  naturally  broader  distribution  tha.n  '‘He.  The 

^  V 

veiccitv  distribution  characteristic  of  this  system  is  —  :::  5%  FWKM 

V 

for  ^He  at  77  ;<. 

E .  Low  Energy  Seam  Production 

The  temperature  of  the  gas  at  the  nozzle  essentially  determines 
the  average  kinetic  energy  of  atom.s  i.n  the  beam.  Production  of  a  very 
Icw-eneriy,  lo:'.g-'.v-ave length  beam  is  therefore  dependent 
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r.C2zle  t:o  very  low  tejiperatures .  This  is  acccT.plished  by  ooolir.g  it 
with  liquid  haliuj:!.  The  problems  of  beam  intensity  and  stability  then 
need  to  be  overcome. 

The  arrangeme-nt  for  cooling  the  nozzle  is  schematically  illustra- 
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ted  in  Figure  7.  It  is  an  adaptation  of  a  "Heli-Tran"  refrigerator 
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and  is  similar  to  the  design  of  Skofronick  and  Pope.  Liquid  helium 
in  the  dewar  is  pressurized  and  a  slow,  steady  stream  forced  to  the 
tip  while  a  portion  of  the  cold  gas  is  fed  back  through  an  outer 
jacket  for  cryogenic  shielding.  The  line  is  also  vacuum-iac.keted  for 
i.mproved  efficiency  and  is  rated  for  0.75  1/hr  average  liquid  heliu.m 
consumption ;  the  cc.nsumption  rate  in  this  application  is  typically 
~  2  1/hr,  probably  due  to  the  high  ther-mai  load  of  the  warm  He  source 
gas . 

The  nozzle  mount  screws  onto  a  copper  "nosepiece"  which  in  turn 
screws  into  the  end  of  the  "cold  finger."  They  are  both  sealed  with 
i.ndium  gaskets  to  promote  good  thermal  contact.  The  source  gas  inlet, 
thermocouple  wires,  a.nd  heater  wire  are  coiled  around  the  cold  fi.nger. 
The  heater  originally  supplied  by  the  manufacturer  was  self-contained 
a.nd  attached  to  the  tip  of  the  cold  finger.  This  did  not  regulate 
well  and  burned  out  frequently.  A  new  heater  wou.nd  directly  onto  the 
nosepiece  is  new  i.nstalled,  a.nd  works  very  well. 

T.he  heater  power  is  regulated  by  a  feedback  circuit  using  the  tip 
thermocouple  signal  to  mai.ntai.n  a  constant  temperature.  The  early 
problems  with  this  method  (probably  poor  thermal  contact  of  the  heater) 
forced  us  to  re.scrt  tc  temperature  control  simply  by  regulating  the 
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flow  of  liquid  helium  until  the  desired  steady-state  temperature  is 
attained.  This  method  worked  fairly  well  but  was  susceptible  to 
long-term  temperature  drifts.  The  most  recent  data  were  taken  using 
the  heater  controller,  which  provides  both  short-term  and  long-term 
temperature  stability  quite  well. 

The  lowest  nozzle  temperature  attainable  with  this  arrangement 
is  =10  K.  A  reasonably  intense  beam  can  be  obtained  only  if  the 
source  pressure  is  less  chan  1  atm.  The  nozzle  also  tends  to  clog 
readily,  probably  due  to  condensation  in  the  aperture.  A  stable  and 
intense  beam  is  produced  at  5=20  K  with  source  pressures  of  -1  y  atm. 
The  velocity  distribution  under  these  conditions  is  =2  j  ^  rWHM, 
only  slightly  poorer  than  the  77  K  beam. 

F.  Electron  Optics 

The  purpose  of  the  electron  optics  in  this  system  is  primarily 
surface  characterization  rather  than  detailed  studies.  It  has  been 
used  in  the  Auger  mode  to  check  the  state  of  surface  cleanliness  and 
in  the  LEED  mode  to  ascertain  the  directions  of  the  surface  lattice 
symmetry  axes.  brief  description  of  it  will  be  given  here. 
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The  electron  optics  head  is  a  PHI  .Model  10-120.  It  consists 
of  an  electron  gun,  four  hemispherical  grids,  and  phosphor-coated 
collector.  It  is  mounted  on  a  bellows  so  that  it  can  be  moved  in 
close  to  t.he  sam.ple  to  wor.k  on  it  and  m.oved  out  of  the  way  when  not 
in  use.  The  electron  beam  is  at  right  a.ngles  to  the  atomac  beam. 
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In  the  LEED  mode,  the  nead  is  controlled  by  a  PHI  Model  11-020 

electronics  system.  The  electron  beam  energy  is  variable  and  the 

hemispherical  grid  system  discriminates  against  inelastically  scattered 

electrons.  In  the  Auger  mode,  the  beam  energy  is  fixed  at  3  keV.  The 
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hemispherical  grids  are  used  as  a  variable  energy  analyzer.  A  PHI 
Model  11-500  Auger  System  Control  sweeps  the  retarding  potential 
through  the  desired  values.  An  ac  modulation  is  superimposed  on  this 
potential  for  lock-in  detection.  Actually,  the  lock-in  monitors  twice 
the  modulation  frequency  in  order  to  take  the  derivative  of  the  elec¬ 
tron  energy  distribution  and  remove  the  background. 
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V .  EXPERIMENTAL  RESULTS 


A.  Alkali  Halides 


The  alkali  halide  (001)  surface  has  been  studied  very  extensively 

g 

by  atomic  beam  scattering.  Boato  et  al.  made  a  thorough  study  of  “He 

diffraction  from  LiF.  Meyers  '  and  Liva"^*^  measured  the  binding 

energies  of  “He  on  LiF  and  NaF,  respectively,  using  the  vacuum  cleaving 

technique  pioneered  by  Houston. Liva  also  experimentally  con- 
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firmed  some  of  the  theoretical  scattering  predictions  of  Cho;.'  and 

Thompson. The  author's  work  on  the  alkali  halides  consists  in  part 

of  an  extension  of  that  work,  employing  many  of  the  same  met.hods  (e.g., 

cleaving  ^  situ) .  More  accurate  binding  energies  were  measured  for 

“He  on  LiF  and  NaF,  and  similar  measurements  were  made  for  ^Ke  on  LiF 
46 

and  NaF,  which  had  not  been  done  before.  This  work  is  summarized 
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in  Section  A.l.  The  band-structure  of  “He  on  NaF  was  verified  and 

matri.x  elements  of  some  higher  Fourier  components  of  the  “He/LiF  atom- 
19 

surface  potential  were  measured.  A  review  of  this  work  is  presented 

in  Section  A,  2.  Data  were  taken  to  test  the  m.ore  recent  scattering 
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theories  for  both  elastic  scattering  and  surface  temperature 
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dependence.  These  studies  are  summarized  in  Section  A. 3. 


The  work  was  done  on  LiF  and  NaF  crystals  cleaved  in  vacuum.  T.he 
surfaces  contam.inate  over  a  period  of  several  days,  and  then  anoti'.er 
cleave  is  made.  Some  of  the  data  were  taken  with  the  surface  cooled 
to  -lOO  K,  but  this  does  not  make  a  dramatic  difference  in  the 
scattering  results. 
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The  surface  structure  of  the  alkali  halide  (001)  plane  is  shown 

in  Figure  8,  along  with  the  conventions  being  used  to  identify  the 

coordinate  axes.  The  rather  corrugated  surface  yields  strong,  sharp 

diffracted  beams.  The  samples  were  bars  about  5  mm  ^  5  mm  x  50  mm 
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obtained  from  .Harshaw  Chemical  Co.  To  promote  cleavability ,  some 
of  the  material  was  hardened  by  y  irradiation  and  some  was  supplied 
doped  with  OH  impurities  ("X-ray  grade"). 


1.  Binding  energies 

As  was  discussed  above,  selective  adsorption  transitions  can  be 
observed  as  anomalies  in  the  specular  beam  intensity.  These  anomalies 
are  generally  minima,  an  example  of  which  is  shown  in  Figure  9.  The 
positions  of  these  features  can  be  used  along  with  Equatior.  (6)  to 
determine  the  binding  energies  E.  of  the  atoms  physiosorbed  on  the 
surfaces.  In  order  to  determine  more  information  about  the  same  inter¬ 
action  potential,  two  isotopes  of  helium  were  used. 

In  the  case  of  LiF,  four  bound-state  energies  were  measured  for 

“He  and  two  for  ^Ke.  Three  energy  levels  were  resolved  for  “he  on 

NaF  and  two  for  ^He.  Several  other  weak  minima  were  observed,  but 

they  could  not  be  resolved  well  enough  to  calculate  energy  values. 
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Preliminary  results  had  been  obtained  previously  '  for  '*He  adsorp¬ 
tion,  as  mentioned  above. 


Data  for  each  energy  level  were  ge.nerally  taken  on  several  cleaved 
surfaces.  Readings  were  also  taken  with  a  variety  of  incident  geo¬ 


metries.  These  are  summariood  i.n  the  form  of  olots  in  K  -K  space 
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for  all  four  systems  studied  in  Figures  10-13.  Most  of  the  data 


He  (l.l  A)  /  LiF 
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were  taken  for  the  1  transition,  with  a  few  ooints  for  G  ,  because 

■5  <  1  ‘  1  »  C 

these  are  the  shortest  reciprocal  lattice  vectors  and  therefore  yield 
the  strongest  selective  adsorption  features.  Care  was  take.n  to  avoid 
collecti.ng  data  for  transitions  which  are  strongly  i.nfluenced  by  band- 
splitting. 

The  average  values  obtained  for  the  bindi.ng  energies  of  ‘‘He/MaF 
a.nd  ^He,  ;.’aF  are  su.trr\arized  i.n  Table  1.  Also  shown  are  the  sta.ndard 
deviations  and  the  nunber  of  detemi.natic.ns  that  were  averaged.  In 
Table  2,  the  sa.-.e  inforr.aticn  for  "He/LiF  and  "He/LiF  is  listed.  The 
sta.ndard  deviations  in  Tables  1  and  2  can  be  compared  to  the  estimated 
■uncertainty  in  the  energy  res'ults.  The  latter  quantity  is  calculated 
as  the  ms  value  of  the  individual  errors  caused  by  'uncertainty  i.n 
the  .T.eas'urer.e.nts  of  d,  ;  and  .k,  .  These  quantities  are  taken  to  be 
dr  =  1;  =  0.2^  and  bk  =  2.23  A“  ‘ ,  resulting  i.n  =  9,03  meV  for  the 

cases  of  interest.  This  is  comparable  to  the  standard  deviations 
found  i.n  the  data.  Use  of  this  data  to  determine  details  of  the  inter¬ 
action  between  helium  atoms  and  alkali  halide  surfaces  is  discussed 
in  Section  VI. A. 


.A  final  word  should  be  said  about  the  cercing  procedures  'used  fo: 
c  and  :.  For  the  gcniomoter  is  ceraed  relative  tc  t.he  symmetry 
directicn.s  of  the  azi.T.uthal  reflection  patterns.  c  is  initiaii'y 
z-erced  relative  to  the  sim.pl-e  holler  face,  but  the  cleaved  crystal 
surface  is  net  nec-cs3ar,.ly  flat  ai.d  parallel  tc  this  face  due  to 
facetting.  This  is  corrensated  for  by  takinc  data  in  four  crvstallo- 
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Table  1. 

Energy  eigen 

values  for  He/NaF 

Isotspe 

3 

E  imeV) 

3 

3td.  dav. 

(T.eV) 

No .  of 

det.ns . 

>4 

he 

0 

-4.92 

0.C5 

7 

'*He 

1 

-1.87 

0.02 

6 

^He 

2 

-0 . 54 

0.05 

3 

^He 

0 

— 1 . 5G 

0.C6 

9 

he 

1 

-1.38 

0.1 

9 

Table  2 . 

Energy  eiger 

lvalues  for  He/LiF 

Isotope 

1 

E.  (f.eV) 

J 

(r.eV) 

No .  of 
det.ns . 

'*He 

0 

-5.90 

0.06 

3 

“Ke 

1 

-2.46 

0.05 

4 

'*He 

2 

-0.73 

0.04 

4 

“■He 

3 

n 

O 

1 

0.02 

4 

^He 

0 

-5.59 

0.08 

5 

^He 

1 

-2.00 

O.Oo 

-* 

...  -.mt,  Sf.'K’*''*'''-  ■> 


43 


The  average  polar  angle  of  the  four  n'.ini.’ta  is  then  used  to  calculate 

c* 


Matrix  elements 


While  the  binding  energies  contribute  information  about  the 
laterally  averaged  aton-surface  potential,  they  tell  nothing  about 
the  periodic  cart,  i.e.,the  G  ^  0  terms  in  Equation  The  periodic 

part  of  the  potential  can  be  probed  by  the  measurement  of  diffracted 
beam  intensities  or  by  determining  matrix  elements  of  the  higher 
Fourier  components.  The  m.any  experi.m.e.ntal  and  theoretical  difficul¬ 
ties  of  the  former  method  make  the  latter  seem  more  attractive.  The 


first  observation  of  such  band-structure  efiec 
and  co-workers  for  Ke  scattering  from  Mar . 


are  tresented  here  for  "He  scattering  from  t.ne 


3  was  .made  by  the  author  ' 

« 

Quantitative  results  1 

(001)  surface  of  LiF .  i 


Figure  14  shews  an  e.xample 
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of  the  sort  of  data  taker,  in  these 
maxi.m.a  in  the  figure  are  selective 
=  0,  G  resonance  i  hereafter 
]  and  the  0(1,1)  resonance.  The  specu- 
function  of  polar  angle  for  several 
that  the  two  maxima  approach  each  other 
decree  apart  in  p.  This  angular  separa- 
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The  0(1,0)  and  0(1,1)  reiaonances  are  strongly  coupled  by  the 

V  (z)  Fourier  component.  The  magnitude  of  the  splitting  can  be  used 
0/1 

to  calculate  the  matrix  element  <0'V  |o>.  The  result  of  this  calcula- 

'  0,1 

19  ^9 

tion,  done  in  collaboration  with  W.  Carlos,  '  is  presented  in  Table 

3.  Also  shown  there  are  the  experimental  matrix  elements  for  all  cases 

which  were  seen  to  give  rise  to  observable  splittings.  Detailed  data 

19 

for  these  matrix  ele.ments  have  been  presented  elsewhere.  These 
matrix  element  values  have  some  interesting  im.plications  for  theoreti¬ 
cal  pair’.vise  S'un  calculations .  This  is  discussed  more  fully  in  Section 
VI,  A. 


3.  Comparison  with  scattering  theories 

It  was  pointed  out  in  Section  III.D  that  the  calculations  of 

,  15  ^  .  ,  16  ,  ...  .  . 

rtears  et  ai.  and  ^arcia  et  a^.,  done  wit.nn  an  eiastic  scattering 

framework,  predict  both  mi.ni.ma  and  .ma-xima  i.n  t.he  specular  intensity. 
.'•Uni.ma  are  predicted  if  the  resonance  state  is  coupled  strongly  to 
both  the  incident  beam  and  an  open  diffraction  channel,  while  ma.xi.ma 
are  expected  for  weak  coupling.  This  was  tested  for  ‘*He  scattering 
from  the  (001)  surface  of  LiF.  Figure  16  is  a  scan  of  specular  inten¬ 
sity  as  a  function  of  azimuthal  angle  at  d  =  70^  (nominal  readi.ng)  . 

The  selective  adsorption  features  seen  in  the  data  are  in  quite  good 
agreement  with  the  theoretical  predictions,  the  strongly  coupled  c0,l) 
states  being  the  deep  minima. 

The  status  of  i.nelastic  scattering  theory  is  on  a  less  firm  basis. 
The  eff-acts  of  the  surface  temperature  in  reducing  elastic  scattering 
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Table  3 . 

‘*Ke,^LiF  matrix  elements  in  meV 

<olv,  Jo> 

0.25  t  0.05 

Jl> 

0.20  r  0.10 

A 

O 

o 

V 

0.10  ±  0.05 

<l|Vi  Jl> 

<0.05 

intensity,  for  example,  is  subject  to  much  debate.  The  standard  methods 

51  52 

employ  a  Debye-Waller  factor  to  account  for  these  losses.  Beeby 

claims  that  the  momentum  transfer  must  be  modified  by  the  acceleration 

53 

due  to  the  attractive  well,  while  Goodman  doubts  that  the  Debye- 
Waller  relation  is  well  established  for  gas-surface  scattering. 

To  examine  these  questions,  a  study  of  "He  scattering  from  the 
(001)  surface  of  WaF  as  a  function  of  surface  temperature  was  under¬ 
taken.  The  major  results  are  stimnarised  in  Figure  17,  which  is  a  plot 
of  the  logarithm  of  the  ratio  of  intensities  at  two  surface  temperatures 
versus  cos'  8.  This  plot  should  be  linear  in  the  context  of  the  Debye- 
Waller  formalism  since 


I(T) 


(10) 


where  I  is  the  intensity  due  to  a  rigid  lattice,  u  is  the  atomic  dis- 
placement,  and  Ik  is  the  momentum  transfer  which  equals  2k  cost)  in 
magnitude  for  the  specular  beam.  It  is  clear  that  such  a  linear  rela¬ 
tion  does  not  hold  for  the  data  in  Figure  17  and  that  a  more  sophisti¬ 
cated  model  is  needed  to  acccu.nt  for  the  effects  of  surface  vibrations 


ntetisity  ratio  for  two  teinperatureis ,  plotted  va . 
test  J-;<juation  (10). 


in  gas-surface  scattering.  If  the  linear  portion  of  the  curve  (roughly 
'3>45^  )  is  interpreted  as  being  a  good  approximation  to  Debve-Waller 
theory,  then  a  surface  Debye  temperature  of  9^  =  411  K  is  calculated. 
From  intensity  versus  surface  temperature  data  at  fixed  angle  (9=45°  )  , 
a  value  9^  =  425  c  20  K  is  obtained.  The  two  values  agree  rather  well, 
and  also  compare  well  with  the  value  of  'Williams  et  al . =  415  K) 
for  He/NaF. 

3 .  Graphite 

Helium  on  the  graphite  basal  plane  is  an  extre.mely  interesting 

system  to  study,  .^s  mentioned  in  Chapter  I,  this  system  has  bee.n  used 

to  carry  out  detailed  thermodynamic  studies  of  physisorbed  films  gaining 

54 

inform.ation  on  phase  transitions  in  two  di.mensions.  Graphite  has  a 
layer  structure  involving  weak  Van  der  Waal's  interlayer  forces.  This 
layer  structure  allows  preparation  of  sam.ples  having  large  surface  area 
(by  means  of  exfoliation)  for  use  in  adsorption  and  phase  transition 
studies.  The  saturated  covalent  bonding  within  the  layers  renders  the 
graphite  basal  phases  chemically  quite  non-reactive,  a  desirable 
property  experim.er.tally  since  it  retards  surface  contamination.  The 

crystallographic  structure  of  the  basal  plane  is  shown  i.n  Figure  15, 
along  with  the  six  lowest-order  reciprocal  lattice  vectors.  T.he 
hexagonal  symmetry  of  t.he  surface  .mesh  is  apparent,  and  the  labelling 
conventions  used  herein  (F.eferer.ce  24  uses  a  dinferent  rcnver.tion) 
are  shown  there.  A  further  convention  is  to  choose  :  =  0°  along  the 

(1,1)  reciprocal  lattice  vector.  Further  motivations  for  employing 

...  ,  ,  .  55. 

grapnite  i.n  scattering  experi.tents  are  its  theoretical  interest  a.nd 

56 

possible  technological  uses,  e.g.  as  a  substrate  for  catalysts. 
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Although  niany  workers  have  used  graphite  for  physisorption  studies, 

employing  a  variety  of  research  techniques , only  Boato  et  al.  besides 

ourselves  have  used  helium  beam  scattering  for  this  purpose.  They  used 

a  63  meV  incident  beam  to  study  the  surface  corrugation  by  means  of 

58 

diffraction  intensities.  They  also  used  a  17.3  meV  beam  to  measure 

59  24 

the  binding  energies  and  matrix  elements  of  '*He  on  graphite. 

The  remainder  of  this  section  presents  the  results  of  the  author's 

study  of  helium  scattering  from  graphite.  In  the  first  subsection, 

methods  of  surface  preparation  and  their  effect  on  the  scattering 

results  are  discussed.  Results  for  the  binding  energies  of  “Ke  and  "He 
26 

on  graphite,  done  at  17.3  m.eV  incident  energy,  are  summarised  in  the 
next  subsection.  In  Subsection  3,  some  interesting  phenomena  associated 
with  level  crossings  are  examined,  at  both  17.3  meV  and  4.7  meV  i.ncident 
beam  energies.  In  the  final  subsection,  experi.mental  evidence  for  in¬ 
elastic  scattering  processes  is  presented,  primarily  at  4.7  meV  incident 
energy  and  lower. 


1.  Sample  preparation 

Samples  for  use  in  atomic  beam  scattering  must  be  natural  single 
60 

crystals.  Grafoil,  which  is  an  exfoliated  form  of  graphite  much  used 
in  thermodynamic  work,  is  not  suitable  for  use  in  beam  scattering  due 
to  its  crystallographic  disorder.  The  same  is  true  even  of  the  pyro¬ 
lytic  graphites,  which  are  highly  ordered  with  respect  to  the  c-a.xis 
orientation  but  completely  disordered  arou.nd  it.  T.he  natural  crystals 
used  in  this  work  are  flakes  about  4  mm  <  4  mm.  ’<  0.1  .mm  in  sits.  They 
do  not  exhibit  a  single  flat  crystal  plane  but  rather  a  mosaic  of 


facets  demarcated  by  fault  lines.  These  facets  are  aligned  well  in 
azimuthal  angle  but  not  in  polar  angle.  The  largest  facets  are  of  order 
1  mm‘  in  area. 

The  crystals  in  their  natural  state  are  embedded  in  a  matrix  of 
calcite  and  silica.  A  quantity  of  this  material,  found  in  Ticonderoga 
(New  York  State)  by  Dr.  T.  S.  Noggle  of  Oak  Ridge  National  Laboratory, 
was  kindly  sent  to  us.  The  graphite  is  extracted  by  attacking  the 
matrix  first  with  hydrochloric  and  then  with  hydrofluoric  acid.  Drs. 

G.  Wagoner  and  .M.  B.  Dowell  of  Union  Carbide  Corporation  and  Dr.  G. 
Felcher  of  Argonne  National  Laboratory  have  also  sent  us  Ticonderoga 
graphite  flakes,  already  extracted,  that  were  used  in  this  work. 

S'' 

The  flakes  are  first  cemented  with  "Vacseal"  “  onto  copper  mounts 
that  fit  into  the  sample  holder.  One  of  two  methods  for  cleaning  the 
surface  are  then  employed.  Most  of  the  data  were  taken  using  surfaces 
which  were  cleaned  by  heat  treatment.  The  sample  is  cleaved  in  room 
air  by  peeling  off  graphite  layers  with  adhesive  tape.  Several  cleaves 
may  be  made  until  a  surface  with  relatively  large  and  well-aligned 
facets  is  obtained.  It  is  then  installed  in  the  scattering  chamber  and 
the  system  is  pumped  down  and  baked.  The  sample  is  heated  to  about 
700  K  for  several  hours  concurrent  with  the  end  of  the  bake.  The  heat¬ 
ing  occasionally  had  to  be  carried  out  twice  before  a  strong  specular 
reflection  was  seen. 

The  second  met.hod  of  surface  preparation  is  cleavi.ng  the  crystal 
in  vacuu.m.  This  is  done  by  affixi.ng  a  thin  cover  glass  to  the  m.cu.nted 
sample  with  epoxy.  After  the  ultra-high  vacuum  is  attained,  t'n.e  glass 


is  detached  along  with  a  layer  of  graphite  by  use  of  a  manipulator. 


The 


scattered  intensities,  linewidths,  and  energy  levels  are  essentially  the 
same  for  both  methods  of  surface  cleaning.  The  scattered  inte.nsities 
are  comparable  to  those  seen  usi.ng  alkali  halides,  but  only  if  the  gra¬ 
phite  surface  is  cooled  to  -100  ;<.  Cooling  gives  rise  to  a  dramatic 
increase  in  elastic  scattering  from  graphite,  as  opposed  to  the  alkali 
halides,  and  all  work  described  herein  was  done  on  a  cold  surface.  T.he 
resonance  linewidths  are  only  slightly  greater  for  graphite  than  for 
alkali  halide  scatteri.ng.  Unlike  the  alkali  halides,  however,  the 
intensities  and  linewidths  were  both  stable  with  time  for  graphite 
using  either  surface  cleaning  technique. 

A  major  disadvantage  of  cleaving  in  situ  is  that  only  one  cleave 
car.  be  made.  Thus,  no  improve.ment  is  possible  (without  opening  the 
chamber)  if  the  facets  are  small  or  badly  misaligned.  The  misalign¬ 
ment  of  the  facets  in  polar  angle  is,  in  fact,  a  proble.m  even  for 
relatively  good  surfaces.  The  magnitude  of  the  misalignment  is 
assessed  by  reflecting  laser  light  from  the  surface.  This  produces 
an  irregular  pattern  of  reflected  light  of  varying  size  dependi.ng 
on  the  quality  of  the  sample.  The  best  samples  produce  a  pattern 
siobtendi.ng  a  solid  angle  of  about  10'“'  sr,  while  the  worst  are  about 
10"^  sr.  For  comparison,  a  polished  copper  sample  mou.tti.ng  plug  gives 
about  10~’  sr.  Thus  it  is  evident  that  even  the  best  crystals  are  far 
from  optically  flat.  The  individual  facets,  however,  are  flat  and  at 
is  only  necessary  that  a  single  facet  be  larger  than  the  beam's 
projection  on  the  surface.  It  is  also  necessary,  of  course,  that  the 
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beara  be  ain-.ed  properly  to  ir-.pinge  on  the  large  facet.  This  is  accom¬ 
plished  by  adjusting  the  height  and  lateral  position  of  the  bean  until 
a  si.ngle,  strong  specular  peak  is  observed.  An  e.v.a.nple  is  shown  i.n 
Figure  19,  obtained  by  scanni.ng  the  detector  through  the  specular  bea.n 
with  fi.xed  i.ncidence  angles. 

Si.nce  the  chose.n  facet  is  not  necessarily  parallel  to  the  sample 
holder,  there  is  still  the  problan  of  zeroing  the  polar  angle.  The 
.nethod  described  in  Section  V.A.l  is  not  feasible  because  the  rotation 
i.nvolved  would  usually  move  the  chosen  facet  out  of  the  bean.  The 
azi.nuthal  zero,  however,  is  detartnined  for  graphite  in  the  sane  manner 
as  for  the  alkali  halides,  described  in  Section  V..^.l. 

The  zero  of  the  polar  angle  is  detem.ir.ed  by  first  writi.ng 


where  c  is  the  true  polar  angle,  i '  is  the  scale  reading,  and  vj  is 
the  scale  readi.nc  whe.n  c  =  O’.  Then,  detector  a.ngle  neasure.nents ,  t_,, 
are  .nude  for  a  set  of  no.ni.nai  polar  angles,  c’.  A  plot  of  versus 
IS  made  and  extranclated  to  t  ,  =  O’.  This  intercept  ef fectivel'.' 


ieternines  the  scai.e  readme 


:rroctuon  v'.  These  plots  turn 


cut  to  be  very  nearly  straignt  li.nes,  as  seer,  in  the  exanple  shew, 
i.n  Figure  20.  Tveueui  xinoar  r’-  rrossicn  ccrrelation  coef  ficie.nts  ’ 


fitter  tron.  u;ittv  cv  or.*/  at 


,  oi  tir.t?  riCwS 


var*/ 


.  1 


0,  a.n  effect  caused  bv  the  various 


misalignments  described  i.n  Toction  IV. 0.  A  detailed  treatm.ent  or  the 
effects  of  apparatus  misalignm.ent  on  the  polar  angle  zeroi;'.g  procedure 
is  given  i.n  the  .Arcendi.x. 
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2 .  Binding  energies 

An  example  of  the  scatter-.ng  data,  demonstrating  selective  adsorp¬ 
tion  minima  in  the  specular  intensity,  is  shown  for  ‘*H&  on  graphite  in 
Figure  21.  Although  the  actual  scattered  intensities  are  not  identical 
for  all  samples  due  to  surface  imperfections,  the  positions  of  the 
minima  are  reproducible  within  experim.ental  error.  This  is  sufficient 
for  present  purposes,  since  t.hese  positions  are  the  input  for  calcula¬ 
tion  of  the  energy  eigenvalues.  The  data  in  this  sectio.n  were  taken 
on  three  surfaces,  one  vacuum-cleaved  and  two  heat-treated,  and  no 
systematic  variations  were  seen  from  sample  to  .sample. 

Five  bound  state  energy  levels  were  resolved  for  “He  on  graphite. 
Energy  levels  were  .measured  for  transitions  via  each  of  the  three 
shortest  reciorocal  lattice  vectors  available,  namely  G  ,  G  ,  and 

0  ,  i  1  I  3 

G  .  The  locations  in  K  -K  space  where  data  were  taken  are  3u.mmarized 
1/1  X  y  * 

in  Figure  22.  Notice  that  the  data  need  be  taken  over  only  a  30° 
interval  in  ?  due  to  the  sixfold  symmetry  of  the  lattice. 

The  average  binding  energies  for  all  “Ke  data  are  summarized  in 

Table  4.  Also  tabulated  there  are  the  standard  deviation  and  n'umfaer 

of  determinations  for  each  energy  level.  The  earlier  results  of 
59 

Boato  et  al .  are  listed  there  for  comparison.  Note  that  there  is  a 
systematic  disagreem.ent  between  the  two  sets  of  eigenvalues.  Although 
this  was  initially  a  source  of  ccncern,  the  discrepancy  car.  be  traced 
to  the  procedure  used  in  zsroi.ng  the  polar  angle  of  the  Genca  appara¬ 
tus ,  and  when  correction  is  made  for  this,  excelle.nt  agreement  with 
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Table  4 . 

Energy 

eigenvalues  for 

** He/  graphite . 

j 

(me'/) 

Std.  dev. 
(meV) 

N'o .  of 
detns . 

Results  of 

a 

Boato  et  al . 
(meV) 

0 

-12. C6 

0.13 

25 

-11.75 

1 

-  6.36 

0.11 

30 

-  6.13 

2 

-  2.35 

0.10 

33 

-  2.65 

3 

-  1.01 

0.09 

32 

-  0.86 

4 

-  0.17 

0.06 

18 

-  -  - 

From  Retarencs  39 


For  ^Ke  on  graphita,  only  chree  energy  levels  were  resolvable. 
This  is  probably  due  to  the  higher  velocity  spread  in  the  ^Ke  incident 
beam.  Its  lighter  mass  causes  Che  ^He  to  be  less  strongly  bound  to 
the  surface,  also  making  it  more  difficult  to  observe  high  j  bound 
states.  The  average  binding  energies  of  the  observed  states  are  pre¬ 
sented  in  Table  5,  along  with  their  standard  deviations  and  the  number 
of  determinations  made  for  each.  The  data  were  taken  over  several 

areas  of  K  -K  soace,  as  shown  in  Ficure  25,  and  were  aaain  taken  for 
X  y  ■  - 

tra.nsitions  induced  by  each  of  the  three  important  reciprocal  lattice 
vectors . 


As  seen  in  Fi 
from  one  region  of 
strongly  coupled  s 
deviations  listed 
uncertainty,  iS .  : 


gures  22  and  23,  no  syste.matic  tre.nds  are  apparent 
K_^-K  ^  space  to  another.  I.n  general,  crossings  of 
tates  were  avoided  in  taking  the  data.  The  standard 
1."  Tables  4  and  5  cc.T.pare  favorably  witii  the  estimated 
0.12  me'/,  calculated  as  in  reoticn  Ti.e  ose 
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Table  5. 

Energy  eigenvalues  for 

^He/graphite 

E  . 

Std.  dev. 

No .  of 

j 

j 

(meV) 

(meV) 

detns . 

0 

-11.62 

0.09 

25 

1 

-  5.38 

0.13 

2  5 

2 

-  1.78 

0.11 

5 

of  t±iese  results  in  detemining  v(z)  for  helium  on  graphite  is  discus¬ 
sed  in  Section  VI. B.  Also  presented  there  is  an  assessment  of  the 
validity  of  the  free-atom  approximiation  used  in  calculating  the  eigen¬ 
values  from  the  resonance  positions  and  Equation  (6) . 


3.  Second-order  resonances  and  splittings 


a.  Using  17.3  meV  incident  beam. 


As  pointed  out  in  Section  III.D/  several  interesting  effects  can 

occur  along  with  band  splitti.ng  when  strongly  coupled  resonances  cross. 

Figure  24  presents  some  specular  intensity  data  plotted  as  a  function 

of  azimuthal  angle  for  various  fi.xed  values  of  polar  angle.  Boatc 
24 

et  al.  observed  a  case  similar  to  this  at  somewhat  higher  energy. 

The  rescna.nce  loci  of  Figure  24  are  clotted  in  K  -K  space  in  Figure 

-  X  y 

25,  along  with  curves  representing  the  free-atom  results.  The  data 

are  shifted  in  polar  ancle  to  compensate  for  facetting  as  in  Section 

V.A.2,  in  this  case,  to  fit  eigenvalues  plus  calculated  band-structure 
29 


corrections . 


Tile  data  in  ricuro  24a  is  clotted  in  region  "a"  of 
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Figure  25,  and  similarly  for  b  and  c.  At  the  lowest  polar  angle,  in 
Figure  24a  and  region  "a",  only  the  pair  of  minima  for  the  0(1,0)  and 
0(0,1)  states  are  visible.  These  are  first-order  transitions,  i.e., 
coupled  to  the  incident  beam  via  first-order  reciprocal  lattice  vectors, 
and  therefore  quite  strong.  As  S  increases,  these  mini.ma  move  inward 
toward  =  0°.  As  this  happens,  another  minimum  appears,  much  weaker, 
at  9  =  72.5°.  It  t.hen  splits  and  moves  outward,  growing  stronger  and 
merging  with  the  other  pair  at  9  =  75.5°.  As  seen  in  Figure  25,  the 
wea.k  minima  correspond  to  the  second-order  transitions  0(1,1)  and 
0(1,1).  The  0(1,1)  state  is  strongly  coupled  to  5(1,0)  by 
and  likewise ,0 (I , 1)  to  0(0,1)  by  Mote  that  the  strong 

coupling  dees  not  occur  between  the  pairs  of  states  that  cross  at 
6  -  75.5°  ;  thus,  there  is  no  obseir/able  splitting  there. 

Consider  next  Figure  24b  and  region  "b"  of  Figure  25.  Here, the 
merged  minima  split  as  9  continues  to  increase,  the  first-order  ones 
movi.ng  inward  to  t.heir  crossing  at  9  ==  77°  and  t.he  second-order  ones 
moving  outward.  Finally,  in  Figure  24c  and  region  "cl’  the  strongly 
interacting  pairs  ceme  together.  In  the  region  of  strong  interaction, 
they  form  two  mixed  states  separated  by  a  band  gap.  After  t.he  crossing 
point  is  passed,  t.he  levels  start  to  approach  the  free-atom  curves 
and  the  second-order  ones  fade  away  as  the  admixture  becomes  weaker. 

This  is  evidently  an  example  of  the  phenomenon  predicted  by 
Weare  et  al.^^  The  second-order  resonances,  not  normally  observable, 
become  manifest  by  "borrowing"  intensity  from  the  first-order  states 
they  are  admixed  with.  From  the  magnitude  of  the  splitting  of  3(1,0) 


and  0(1,1),  the  matrix  element  <OiV  0>  is  calculated  to  be  0.25cC.03 

0  t  1 

24 

meV.  These  results  are  consistent  with  these  of  Boa to  et  al.  (see 
also  Taible  11)  . 

b .  Using  4.7  meV  incident  beam. 

The  use  of  a  very  low-energy,  long-wavelength  incident  beam  has 

several  desirable  features.  The  constraints  imposed  by  the  kinematics 

of  the  scattering  result  in  fewer  open  diffraction  channels  at  these 

energies.  The  number  of  selective  adsorptions  allowed  also  decreases 

a.nd  the  resonances  spread  cut  in  space.  Thus,  the  situation  is 

somewhat  simplified  as  compared  to  the  17.3  .meV  beam.  An  application 

of  this  low  energy  beam  is  presented  here  for  helium  scattering  frem 

graphite,  namely  a  study  of  some  splittings,  which  are  also  slightly 

more  spread  out  in  K  -K  soace. 

X  y  * 

A  situation  similar  to  that  just  previously  described  is  seen 

in  Figures  26  and  27,  t.his  data  taken  with  4.7  meV  incident  energy 

0—1 

corresponding  to  -  3,0  A  .  The  region  being  examined  again  in- 
eludes  both  first-  and  second-order  transitions,  the  first-orders  being 
j  =  1  states  this  time.  For  the  lowest  values  of  9,  the  minimum  at 
the  higher  corresponds  to  the  1(0,1)  state,  while  the  minimum  at 
lower  0  corresponds  to  0(1,1).  As  9  increases,  both  resonances  move 
inward  toward  their  .nominal  crossing  poi.nts  at  1  =  0^.  .At  t  -  SO'’, 
the  0(1,1)  beccm.es  degenerate  at  5  =  0°  wit.h  t.he  0(1,1)  symmetrical 
to  it  (.not  shown  in  Figure  26  for  simplicity).  The  1(0,1)  resonance, 
t  off  by  its  admixture  with  the  0(1,1)  and  dees 
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situations  of  high  symmetry,  but  there  is  no  obvious  symmetry  here. 

Thera  is  some  evidence  that  it  may  be  caused  by  inelastic  effects. 

At  any  rate,  an  estim.ate  of  the  matrix  element  for  the  splitting  can 

be  calculated  by  extrapolating  the  observed  data.  It  is  found  to  be 

<llv  |o>  =  0.16  3  0.03  meV,  which  again  agrees  within  excerimental 
0/1 

24 

error  with  the  result  of  Boato  et  al.  obtained  at  the  higher  beam 
energy. 

The  (1,0)  states  are  only  weakly  coupled  to  the  (1,0)  and  (1,1) 
states  they  cross.  For  j  =  4  and  j  =  2, they  simply  cross  with  no 
noticeable  disturbance,  as  expected.  The  3(1,0)  transition,  however, 
cresses  through  the  1 (1,0) -0 (1,1)  splitting  and  it  also  disappears 
in  this  vicinity.  Since  the  coupling  in  this  case  is  weak,  the  situa¬ 
tion  is  somewhat  different  than  for  the  (1,0)-(1,1)  mixture.  Neither 
case  is  really  well  understood  theoretically.  Finally,  notice  that 
the  average  binding  energies  of  the  (1,0)  resonances  i.n  this  region  are 
-O.lrneV  higher  than  the  overall  averages  previously  established.  This 
is  seen  to  be  the  case  with  the  17.3  meV  beam  energy  also. 


4.  Inelastic  effects 

Even  at  17.3  meV  i.ncident  energy,  there  is  some  evidence  of  the 
presence  of  inelastic  processes  operating  in  the  scattering.  The 
experimentally  observed  selective  adsorption  signatures  are  virtually 
all  intensity  minima,  whereas  the  elastic  scattering  theories  of  Wearv 
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et  al.  and  darcia  et  al .  predict  maxima  in  a  few  case.s.  Mcr', 
lin;;  evidenje  ca.n  be  adduvjed  bv  use  of  the  lower  enerqv  beam. 
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Figure  30  shows  specular  intensity  data  as  a  function  of  polar 
angle  at  ■-{)  =  0®  for  two  different  incident  energies.  Figure  30a  is 
data  for  the  beam  energy  at  16.8  meV,  while  Figure  30b  is  for  the 
5.1  meV  beam  energy.  The  data  were  taken  on  the  same  surface  and 
incidence  conditions  other  than  experimental  variables  were  held  as 
constant  as  possible. 

The  major  qualitative  trends  in  the  data  can  be  explained  withi.n  an 
elastic  scattering  framework.  In  Figure  30a  the  intensity  in  the  specular 
beam  shows  an  overall  decreasing  trend  with  6.  This  is  probably  caused 
by  the  intensity  goi.ng  into  the  diffraction  channels  which  become  open  at 
lower  9.  When  the  beam  energy  is  5.1  meV,  fewer  diffraction  channels 
open  and  this  is  reflected  in  the  relatively  flat  (except  resonances) 
intensity  scan  of  Figure  30b.  Si.nce  there  are  fewer  open  cha.nnels, 
however,  it  is  expected  that  the  specular  intensity  at  5.1  meV  be 
greater  than  that  at  16.3  meV  for  any  polar  angle.  Prom  Figure  30  it 
is  seen  that  the  intensities  are  about  the  same  magnitude.  This  fact 
also  contradicts  t.he  predictions  of  a  Debye-Waller  model.  A  highly 
probable  explanation  of  this  is  the  presence  of  inelastic  scattering 
processes  draining  inte.nsity  from  the  specular  beam. 

In  order  to  perform  a  more  direct  test  of  this  hypothesis,  the 
beam  energy  is  further  lowered  to  about  2.9  meV.  At  this  incide.nt 
energy,  only  one  diffracted  beam,  the  (1,1)  is  allowed  at  =  0°, 
and  it  is  an  in-plane  beam.  Thus,  it  is  possible  to  monitor 
the  total  diffracted  flux  from  a  well-defined  arystal  plane  a.;  the 
polar  ar.cl;  is  scanne-d  through  a  selective  adsorption  resonance. 
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Any  unitary  elastic  theory  predicts  that  an  intensity  minimuin  in  one 
of  the  beams  should  give  rise  to  a  maximum  in  the  other.  Results  for 
the  1 C1,Q) -1  (0 , 1)  transition  are  presented  in  Figure  31,  where  it  is 
seen  that  minima  in  both  beams  occur.  It  is  felt  that  this  is  an 
unambiguous  demonstration  of  the  importance  of  inelastic  scatteri.ng, 
and  probably  as  well  as  can  be  done  without  energy  analysis  of  the 
scattered  beams. 


Diffraction 
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VI .  DISCUSSION 

A.  Alkali  Halides 

1.  Laterally  averaged  potentials 

A  major  experimental  result  of  Section  V.A  is  the  eigenvalue 

spectrum  of  helium  adsorbed  on  NaF  and  LiF  surfaces.  These  values  are 

the  input  for  the  generation  of  semi-empirical  potentials  V  (z) .  Two 

0 

isotopes  of  helium  are  used  to  increase  the  input  available  for  these 

calculations.  The  eigenvalue  data  also  would  serve  as  a  test  of  any 

first-principles  calculation  of  V  (z) ,  but  this  has  not  yet  been  done. 

0 

Several  semi-empirical  model  potentials  have  been  used  with  success 

by  Cole  et  al,^*^  Two  of  these,  the  3-9  potential^*^ '  and  the  shifted 

64 

Morse  hybrid  potential,  have  been  fit  to  the  alkali  halide  eigen¬ 
value  data  and  are  discussed  here.  The  3-9  ootential  can  be  written 


where  is  the  well  dept.h,and  j  is  the  distance  for  which  the  potential 
is  equal  to  zero  (Zj  simply  shifts  t.he  z-a.xis  without  changing  the 
eigenvalue  spectrum) .  The  values  of  the  parameters  are  summarized  in 
Table  6,  and  the  calculated  eigenvalues  are  compared  to  experi.ment  in 

Table  7.  Equation  (12)  has  the  correct  asymptotic  form  of  — ,  far  from 

66  1 
the  surface,  and  the  coefficient  C,  of  the  asvmototic  —  can  be  cal- 

3  ‘  _  3 

culated  from  D  and  J.  The  C,  calculated  in  this  wa'.- ,  however,  is 
e  3  .  '  ' 

about  twice  the  theoretical  value  of  Bruch  and  Watanabo .  This 
discrepancy  is  probably  due  to  the  fact  that  the  data  is  representative 
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Tatile  6A.  Parameters  of  the  3-9  potential  (from  Ref.  64). 


c 

C3 

System 

(meV) 

(i) 

3  3 

(meV 

(A) 

‘D 

He/NaF 

7.64 

1.91 

138.3 

0.59 

3.58 

He/LiF 

8.92 

1.94 

169.2 

0.70 

3.93 

Table  6B.  Parameters  of 

the  shifted  Morse 

hybrid  potential 

( from 

Ref . 

64)  . 

D 

System  (m.eV) 

(A  ) 

1 

z 

(A) 

D(1  +  1) 
(A)  (meV) 

a 

C. 

meV-A- 

Ke/NaF  6.5176 

1.2269 

0.05937 

4.263 

1.699  6.905 

72 . 5 

He/LiF  7.740 

1.2929 

0.075 

3.432 

1.634  3.321 

SI.  7 

^Theoretical  value,  from  Reference  67 
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Table  7. 

Comparison 

of  experimental 

and  calculated 

eigenvalues.  )  \ 

System 

1 

1 

Ej  1  (expt. ) 

(meV) 

|Ej 1  (3-9)^ 

]  ' 

(meV) 

(meV) 

^He/NaF 

0 

4.50 

4.53 

4 

4.558 

1 

1.38 

1.40 

1.380  K 

2 

— 

0.33 

0.245  : 

3 

— 

0.05 

V 

0.013  p 

f  i 

“He/NaF 

0 

4.92 

4.89 

k 

4.338 

1 

1.37 

1.82 

1.883  1'? 

2 

0.54 

0 . 56 

0.471  ^3 

3 

— 

0.13 

0.071 

^He/LiF 

0 

5.59 

5.56 

n 

5.616  IJ 

1 

2.00 

1.94 

1.945  |j 

2 

— 

0.54 

0.412  ^ 

3 

-- 

0.11 

0 . 045 

t* 

“He/LiF 

0 

5.90 

5.94 

5.940 

1 

2.46 

2.44 

2.437  !' 

2 

0.78 

0 . 36 

4 

0.731  I*' 

3 

0.21 

0.24 

0.141 

4 

— 

0.C5 

0.014 

1 

^From  Reference  64 


of  acoms  localized  near  the  surface  rather  than  in  the  asymptotic 

region.  A  convenient  property  of  Equation  (12)  is  that  it  predicts, to 

6>  5 

very  good  approxi.mation ,  an  analytic  eigenvalue  spectrum; 


n  =  n.  -  a  1e, 


where  1  is  a  mass  reduced  quantum  number  allowing  utilization  of  data 

1  !/'> 

for  both  isotopic  forms  and  is  defined  by  n  =  ( j  r  /m 

The  shifted  Morse  hybrid  potential  is  based  on  the  idea  that  the 
Morse  potential  is  a  good  approximation  to  V,^  (z)  near  t.he  minimum  but 
has  the  incorrect  asymptotic  form.  Thus,  for  z  greater  than  some  poi.nt 

Cj 

z  ,  the  form  —  is  assumed  with  C  bei.ng  the  theoretical  value.  For 

^  7. 

z  less  than  z  ,  a  shifted  Morse  function 
P 

=  D(x'  -  2x  -  d),  X  =  (14) 

is  used,  where  D  is  the  well  depth,  d  is  the  shift,  and  z^  is  the 
equilibri'um  position.  The  parameters  of  equation  (14)  are  fitted  to 


the  energy  eigenvalues  and  the  functions  are  matched  at  z  .  Tables 

P 

6  a.nd  7  give  t.he  param.citers  and  predicted  eigenvalues  of  t.he  shifted 
Mores  hybrid  potential  for  helium  on  the  alkali  halides.  .4  plot  of 


t.hese  potentials  is  give.n  in  Figure  32  for  helium  on  NaF. 


2.  Surface  periodicit 


The  matrix  elements  m.easured  for  He  on  biF  can  be  com.parec.  to 
soretical  calculations  done  by  T3uchida°°  and  by  Chicw  and  Th.omcscn. 
3  method  they  emtlov  is  to  model  the  cote.ntial  as  a  su.d  of  -air 


me  thod 
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Table 

8. 

Experimental  and  theoretical 
elements  (in  meV) . 

values  of 

“He/LiF  matrix 

Chow  and 

Thompson 

b 

Tsuchrda 

Experimental 

io> 

0,32 

0,43 

0.25  ±  0.05 

|1> 

0.14 

0,20 

0.20  ±  0.10 

A 

O 

*-> 

|o> 

0.21 

0.28 

0.10  ±  0.05 

A 

< 

!i> 

0.09 

0.14 

<0.05 

^From  Reference  17 
b 

From  Reference  63 
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interactions  between  the  He  atom  and  the  ions  of  the  crystal,  and  they 
compute  the  Fourier  expansions  of  the  resulting  expressions.  From  the 
one-dimensional  Schroedinger  equation  using  (z) ,  the  wave  functions 
are  calculated.  The  theoretical  matrix  elements  are  then  found  by 
numerical  integration.  Tsuchida  represented  the  potential  as  a  sum  of 
12-6  pair  interactions,  while  Chow  and  Thompson  used  sums  of  Yukawa  -6 
pair  potentials.  The  results  of  their  calculations,  along  with  the 
experimental  values  for  comparison,  are  presented  in  Table  8. 

Notice  that  the  theoretical  models  overestimate  the  strength  of 

29 

Vj  j  as  compared  with  ^ .  Carlos  has  pointed  out  that  this  may  be 
due  in  part  to  the  approximations  made  in  the  pair  potentials.  In 
particular,  the  He-Li  pair  interaction  is  approximated  by  that  for 
He-He,  and  He-F  interaction  is  approximated  by  that  for  He-Ne.  This 
is  consistent  with  the  He/LiF  diffraction  data,  which  can  be  inter¬ 
preted  quite  well  using  only  the  (0,1)  component  of  a  surface  corruga¬ 
tion  function. 


B .  Graphite 

1.  Application  to  physisorption 

As  in  the  case  of  the  alkali  halides,  the  eigenvalues  constitute 

29  69  70 

the  spectroscopic  data  of  the  adatom  states.  Carlos  and  Cole  '  ' 

have  fit  several  model  potentials  to  data  for  both  isotopes.  These 

include  pairwise  sums  of  12-6  potentials,  12-8-6  potentials,  exponential 

-  6  potentials,  and  Yukawa  -  6  potentials.  In  all  cases, they  find  that 

the  Ej  can  be  fit  to  these  potentials,  but  the  matri.x  elements  measured 

24 

for  He  on  graphite  by  Boato  et  al.  ca.nnot  be  accounted  for.  This 
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situation  is  remedied  by  the  introduction  of  anisotropy  into  the  poten¬ 
tial  corresponding  to  the  anisotropy  of  the  bonding  orbitals  in  the 
graphite  basal  plane.  They  sum  over  pair  potentials  of  the  form 


U(r) 


A  e 


-ar 


(1 


+  cos 


6)- 


3a  dC, 


TTr° 


f 


cos" 


(15) 


and  obtain  a  laterally  averaged  potential 

r  n  CO 


V,  (z)  » 

’  aa 


-az 

e  -t-Y.azEj  (az) 


-SCjd  (z+nd) 
n=0 


(16) 


where  a^  is  the  unit  cell  area,  d  is  the  interplanar  spacing,  and 

is  an  exponential  integral.  The  parameters  of  (z)  are  fit  to 
the  eigenvalues  tabulated  in  Section  V.B.  The  pairwise  sums  also  yield 
t.he  higher-order  V-(z),  which  are  fit  to  the  matrix  elements.  The  para- 
meters  of  the  modified  Yukawa  -  6  potential  are  listed  in  Table  9,  while 
the  eigenvalues  and  matrix  elements  calculated  from  it  are  compared  to 
experimental  values  in  Tadjles  10  and  11.  The  agreement  is  very  good, 
and  this  is  thought  to  be  a  good  approximation  to  the  helium-graphite 
interaction.  On  the  other  hand,  the  optimum  value  of  C^,  given  in 
Table  9,  is  about  half  the  theoretical  value. This  is  probably  be¬ 
cause  the  atoms  in  resonance  do  not  sample  much  of  the  asymptotic 
region.  A  plot  of  the  potential  over  various  lattice  sites  is  shown 
in  Figure  33. 

A  further  test  of  both  experiment  and  theory  is  to  compare  these 
results  with  the  data  from  experiments  on  thermodynamic  properties  of 
He  films  adsorbed  on  graphite.  Before  this  can  be  done,  however,  the 
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Table  9.  Parameters  of  modified  Yukawa-6  potential,  from  Reference  69 


Y 


R 


a 


c 

3 


A 


0.4  -0.29  3.25  i"*  83  meV-i^  47.2  eV 


a 

Fixed  a  priori  from  dielectric  properties  of  graphite 


Table  10. 

Experimental 

He/graphite 

and  theoretical 
(in  meV) . 

energy  eigenvalues 

of 

n 

‘‘He/graphite 

Experimental  Theoretical 

^He/graphite 

Experimental  Theoretical 

0 

-12.06 

-12.12 

-11.62 

-11.61 

1 

-  6.36 

-  6.36 

-  5.36 

-  5.34 

2 

-  2.85 

-  2.77 

-  1.78 

-  1.88 

3 

-  1.01 

-  0.94 

— 

-  0.44 

4 

-  0.17 

-  0.22 

— 

— 

a 

From  Reference  69 


i 
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Table  11. 

Theoretical  and  experimental 
ite  (in  meV)  . 

matrix  elements  for  He/graph- 

Theoretical^ 

Experimental 

Boato  et  al.^  Penn  State 

<o|v  lo> 

0  1 

-0.27 

-0.28  -0.25 

<i|v  |l> 

0  1 

-0.21 

-0.185 

<2iv  |2> 

0  1 

-0.13 

-0.12 

<3jv  |3> 

0  1 

-0.06 

-0.08 

A 

o 

> 

O 

V 

0.20 

0.195  0.19*^ 

<o|v  |2> 

0  1 

0.13 

0.125 

<o|v  |3> 

0  1 

0.08 

0.09 

<o!v  |4> 

0  1 

0.04 

0.03 

<l|v  \2> 

0  1 

0.16 

0.16 

<1|V  |2> 

0  1 

0.10 

0.10 

<2|v  1 3> 

0.09 

0.11 

0  1 


a 

From  Reference  69 

b 

From  Reference  24 
c 

Average  value 


free-atotn  binding  energies  must  be  modified  to  account  for  band- 

structure.  A  full  band-structure  calculation  was  done  by  Carlos  and 

Cole,  '  using  the  data  of  the  Genoa  group  and  of  Section  V.B. 

The  result  of  their  calculation  at  scattering  energies  demonstrates 

the  validity  of  the  free-atom  approximation  used  to  obtain  the 

(far  from  level  crossings) .  Their  calculation  in  the  energy  regime 

of  physisorbed  sub-monolayer  He  films  shows  a  greater  degree  of  corruga- 

72 

tion  than  was  previously  believed,  implying  somewhat  less  free-particle 
character  for  the  adatoms.  This  result  reflects  the  anisotropy  of  the 
potential . 

The  computed  band-structure  corrections  at  ic  =  0  are  -  0.16  meV 
for  *He,3nd  -  0.11  meV  for  ’He.  This  yields  E^  =  -12.22  neV  for  ‘‘He, 
and  Ej=-ii.73  .meV  for  ’.He  as  t.he  ground  state  energy  of  a  single  atom  ad¬ 
sorbed  on  graphite.  The  ground  state  energy  per  atom  in  an  adsorbed 
film  is  equal  to  the  chemical  potential  Uj  at  T  =  OK.  This  quantity 
has  been  derived  from  thermodynamic  data  by  Greif  et  al.^^  for  ’He 
and  ‘‘He  films  on  graphite.  After  correction  is  made  for  the  latent 
heat  of  “He^"*  (  ’He  is  not  e.xpected  to  condense^^)  ,  they  obtain 
Uj,  •  -12.27  t  0.2  meV  for  “He,  and  =  -11.72  t  0.2  meV  for  ’He.  These 
values  are  in  excellent  agreement,  which  seems  significant  considering 
the  vast  differences  between  the  two  types  of  experi.ment . 

2.  Evaluation  of  scattering  theory 

It  is  clear  from  the  previous  treatment  that  the  kinematical 
analysis  has  a  firm  fo'undation.  Results  of  dynamical  scattering 
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theory,  however,  are  still  subject  to  doubt.  Considering  the  success 
of  elastic  theories  in  predicting  resonance  lineshapes  for  He  scattering 
from  LiF  and  their  partial  success  for  graphite,  it  seems  that  inelastic 
processes  may  be  the  cause  of  discrepancies. 

A  major  success  of  elastic  theory  is  the  prediction  by  Weare 
et  al.^^  of  intensity  borrowing  in  level  crossing.  The  data  of 
Section  V.B  demonstrate  that  this  prediction  is  quite  consistent  with 
experiment  for  both  moderate  and  very  low  beam  energies.  The  same 
theory,  however,  predicts  selective  adsorption  intensity  maxima  which 
are  observed  e.xperiinentally  to  be  minima.  The  conjecture  that  this  is 
caused  by  inelastic  scattering  is  consistent  with  the  large  increase 
in  signal  produced  by  lowering  the  surface  temperature.  This  is 
intuitively  sensible  since  weak  interlayer  forces  in  graphite  can  be 
expected  to  produce  soft  vibrational  modes. 

The  data  also  show  a  resonance  disappearing  in  the  vicinity  of  a 

level  crossing,  a  phenomenon  that  occurs  at  both  moderate  and  low 

beam  energies.  This  is  also  unexplained  by  otherwise  successful 

6  3 

theories.  Weare  has  reproduced  this  theoretically  by  using  an 
optical  potential  to  incorporate  inelastic  effects  into  the  scattering 
formalism.  Chow^^  has  also  used  the  optical  potential  to  show  that 
elastic  intensity  maxima  are  transformed  into  minima.  The  optical 
potentials,  unfortunately,  provide  virtually  no  detail  about  the 
atom-surface  interaction. 

Finally,  the  data  taken  with  a  very  low  energy  incident  beam  show 
clear  and  direct  violations  of  any  unitary  elastic  theory.  Thus, 
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despite  its  valuable  successes,  atomic  beam  scattering  theory  must 
directly  confront  the  problem  of  inelastic  scattering  mechanisms  in 
order  to  progress  much  farther.  At  the  very  least,  it  should  be  able 
to  predict  the  conditions  for  which  an  elastic  framework  is  appropriate. 

C.  Conclusions 

The  major  results  of  this  wor.k  are  the  measurement  of  binding 
energies  of  ^He  and  “He  on  NaiF,  LiF  and  graphite,  the  measurement  of 
several  matri.x  ele.me.nts  of  “He  on  LiF  and  grapnite,  and  t.he  examina¬ 
tion  of  resonance  lineshapes  for  “He  scattering  from  LiF  and  graphite, 
t.he  latter  at  several  incident  beam  energies.  These  results  have 
several  interesting  implications  for  both  atomic  beam  scattering  theory 
and  t.he  physics  of  two-dimensional  systems.  In  general,  it  seems  that 
an  elastic  framework  is  more  appropriate  for  He  scattering  from  alkali 
halide  surfaces  t.han  from  the  basal  plane  of  graphite.  Scattering 
studies  on  graphite,  especially  at  very*  low  incident  energies,  demon¬ 
strate  the  presence  of  important  inelastic  processes.  Other  experi¬ 
ments,  however,  e.xhibit  phenom.ena  at  both  moderate  and  very  low  beam 
energies  that  are  consistent  with  an  elastic  formalism,  such  as 
intensity  borrowing.  The  problem  of  determining  the  conditions  under 
which  an  elastic  framework  is  appropriate  remains  to  be  solved.  The 
data  presented  here  indicate  the  usefulness  of  the  kinematical  approach 
regardless  of  the  status  of  dynamical  theory.  Within  a  kinematical 
framework,  moreover,  the  free-atom  approximation  is  seen  to  be  good 
far  from  level  crossings. 
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The  free-atom  approximation  is  used  herein  to  calculate  the  bind¬ 
ing  energies  of  ‘‘He  and  ^He  on  Naf,  LiF,  and  graphite.  From  these,  a 
set  of  laterally  averaged  potentials  (z)  are  generated  for  these 
systems.  (z)  is  well-oharacterized  for  heli'om  on  both  alkali  halides 
and  graphite.  There  are  still  problems,  however,  with  these  models. 

In  particular,  the  coefficient  of  the  potential^  asymptotic  tail  differs 
greatly  from  theory  in  all  these  cases. 

.Matrix  elements  of  the  V-(z)  were  also  measured  for  “He  on  LiF 

G 

at  17.3  meV  and  '*He  on  graphite  at  4.7  meV.  The  LiF  measurements 

indicate  that  the  ^  component  dominates  the  surface  periodicity. 

In  the  case  of  graphite,  the  matrix  elements  measured  at  4,7  meV 

24 

beam  energy  agree  within  error  with  those  of  Boato  et  al,  measured 

69 

at  a  higher  energy.  The  latter  were  employed  by  Carlos  and  Cole, 
along  with  the  “He  and  ^He  eigenvalues,  to  generate  the  form  of  the 
heli'om-graphite  potential.  It  is  found  that  an  anisotropic  form  is 
required.  They  then  use  the  potential  to  calculate  the  band-structure 
of  helium  on  graphite. 

The  major  conclusion  of  these  calculations,  based  partly  on  the 
data  presented  herein,  is  that  the  adsorbed  He  atoms  feel  a  larger  corru¬ 
gation  than  was  previously  thought  to  exist.  This  has  obvious  implications 
for  the  thermodynamics  of  helium  films  on  graphite.  For  the  comparison 
which  has  been  made  already,  the  scattering  data  and  thermodynamic 
data  are  in  remarkable  agreement  for  the  ground  state  energy  of  both 
“He  and  ^He  on  graphite.  Efforts  are  continuing"^*^  in  the  calculation 


of  thermodynamic  properties  from  scattering  results  on  this  system. 
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D.  Suggestions  for  Future  Work 

In  view  of  the  isiportant  role  played  by  inelastic  scattering 
aiechcinisms ,  a  high  priority  should  be  given  to  energy  analysis  of  the 
scattered  beams.  This  can  be  done  by  the  time-of-f light  technique. 

It  will  provide  importa.nt  information  about  the  surface  vibrational 
inodes,  and  will  probably  contribute  to  the  u.nderstanding  of  the 
scattering  mechanisms  as  well. 

Experiments  should  be  done  on  a  variety  of  other  surfaces.  Semi¬ 
conductors,  metals,  diamond,  and  secondary  cleavage  planes  of  the 
alkali  halides  are  a  few  of  the  possibilities.  Another  interesting 
and  valuable  experiment  would  be  scattering  of  helium  from  adsorbed 
overlayers  of  the  heavy  rare  gases.  Rars-gas-plated  graphite  basal 
planes  are  a  particularly  interesting  possibility. 
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APPENDIX 

Misalignment  and  Angular  Measurement 

Al.  Introduction 

The  determination  of  the  incidence  angles,  9  and  'p,  is  a  difficult 
procedure  if  high  accuracy  is  required.  One  problem  is  caused  by  the 
crystal  surface  being  composed  of  facets  misaligned  in  polar  angle. 
Another  problem  is  caused  by  the  apparatus  misalignments  outlined  in 
Section  IV. C  ;  namely, the  beam,  the  surface,  and  ail  rotation  axes  not 
being  properly  positioned  with  respect  to  each  other.  The  two  problems 
are  not  independent  of  each  other, and  the  occurrence  of  both  together 
is  extremely  difficult  to  treat.  Because  the  facetting  is  not  very  bad 
in  the  alkali  halides,  the  averaging  procedure  described  in  Section  V.A.l 
is  adequate  if  care  is  taken  to  minimize  the  degree  of  misalignment. 

The  relative  flatness  of  the  surface,  in  turn,  makes  the  alignment  pro¬ 
cedure  easier  to  perform.  Work  with  the  graphite  surface,  however, 
entails  a  more  severe  set  of  problems. 

A2 .  Angle  .Measurements  for  Graphite 

The  facets  of  the  graphite  surface  typically  are  small  and  badly 
misaligned  (see  Section  V.3.1).  It  is  clear  that  the  averaging  techni¬ 
que,  appropriate  for  the  alkali  halides,  is  not  suitable  for  determining 
0  in  graphite  work  since  an  entirely  new  facet  will  probably  be  sampled 
upon  rotation  by  60“  to  a  crystallographically  equivalent  position. 

Instead,  the  procedure  for  determining  0  described  in  Section  V.3.1 
is  employed.  This  assumes  that  the  polar  angle  scale  reading  0'  is 
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equal  to  the  true  polar  anqle  6  except  for  a  constant  correction  vihich 

can  be  found  by  zeroing  the  angle  scale.  Detector  angle  measurements, 

6^  ,  are  made  for  a  set  of  nominal  polar  angles  0'  and  are  extrapolated 

to  0,  =  0.  This  yields  the  value  of  0'  at  0  =  0®  which  can  then  be 
d 

used  to  correct  all  9'  measurements.  This  procedure,  however,  is  also 
affected  by  the  misalignments  and  its  accuracy  is  therefore  not  neces¬ 
sarily  assured.  Thus,  an  assessment  of  the  effects  of  misalignment  on 
the  polar  angle  zeroing  procedure  was  carried  out. 

A3 .  Effects  of  Misalignment 
A3a.  Surface  Mormal  Misalignments 

If  the  surface  normal  is  not  collinear  with  the  azimuthal  axis 
(a  tilted  facet,  for  example) ,  the  polar  angle  changes  as  the  azimuth 
is  varied.  Therefore, data  are  taken  only  within  a  few  degrees  of  the 
({)  value  at  which  the  zeroing  procedure  is  done.  A  new  zeroing  run  is 
performed  when  data  are  wanted  at  a  different  >{).  The  tilt  will  also 
introduce  some  error  into  the  detector  angle  reading,  but  this  error 
enters  linearly  to  good  approximation  since  the  tilt  angle  is  small  and 
therefore, it  is  taken  care  of  by  the  extrapolation  procedure. 

A3b.  Rotation  Axes  Misalignments 

The  misalignments  of  the  polar  and  detector  angle  rotation  axes 
are  difficult  to  treat  analytically,  so  an  empirical  approach  was 
taken.  In  one  case,  for  example,  a  zeroing  run,  0'  vs.  0^  ,  was  first 
carried  out.  Then, the  detector  angle  axis  was  moved  =3  mm  laterally. 
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The  same  zeroing  run  was  then  done  again  resulting  in  different  detector 
angle  readings  but  having  essentially  the  same  slope  and  the  same  value 
of  These  results  are  presented  in  Table  12.  Although  it  is  not 

possible  to  draw  general  conclusions,  no  inconsistencies  were  found  to 
result  from  these  misalignments  and  it  is  at  least  reasonable  to  expect 
the  extrapolation  procedure  to  compensate  for  their  effects. 

A3c.  Beam  and  Surface  Misalignments 

The  major  remaining  problems  are  the  baam  and  the  sample  surface 
not  being  on  the  polar  angle  rotation  axis.  Only  these  two  misalign¬ 
ments  are  considered  here  to  simplify  the  analysis.  First,  define  d 

B 

as  the  distance  between  the  beam  and  the  axis  of  rotation  d^  as  the 
distance  between  the  surface  and  the  axis,  and  R  as  the  radius  of  the 
detector  swing.  Let  a  =  90 *-0,  the  angle  between  the  beam  and  the  sur¬ 
face,  so  that  in  the  ideal  case  of  no  misalignment  '5^  =  2a.  In  the  case 
of  misalignment  being  present,  it  is  found  by  geometrical  construction 
that 


(Al) 


This  transcendental  equation  is  solved  numerically  and  it  is  found  that 
the  extrapolation  procedure  is  valid,  but  only  for  the  small  degree  of 
misalignment  maintained  throughout  the  experiment. 


R  Sin  y  ,  -  d 


tan  2a  = 


B 


R  cos  0  ,  - 
d 


"B 


tana  sina 
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Table  12.  Zeroing  results  before  and  after  detector  axis  change. 


Data  of 

4>  * 

9/12/79 

17<* 

Data  of  9/15/79 
■+>  =  17“ 

6' 

9 

d 

9' 

®d 

direct  beam 

0.54 

direct  beam 

0.80 

80 

20.78 

80 

21.06 

75 

30.83 

75 

30.58 

70 

40.53 

70 

40.74 

65 

50.54 

65 

50.75 

60 

60.25 

60 

60.33 

55 

70.01 

55 

70.00 

50 

79.93 

50 

79,72 

45 

89.54 

45 

89.43 

40 

99.32 

40 

99.22 

35 

109.11 

35 

109.08 

30 

118.85 

slope  =  -1.96  slope  =  -1.95 

9'  =  0.42-  0'  = 

“  0 


0.39 
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